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SUMMARY 


The  transient-response  calculation  method  reported  in 
USAAVLABS  Technical  Report  65-65  ia  extended  to  yield  aerodynamic 
loadinga  and  elastic  -blade  motions  in  response  to  inputs  of  shaft  torque. 
Whereas  transient  blade -load  calculations  were  previously  restricted  to 
low-advance-ratio  flight  conditions  in  which  the  flow  relative  to  the  blade 
is  primarily  two-dimensional,  aerodynamic  loadings  are  also  calculated 
in  this  study  for  flight  conditions  in  which  the  flow  relative  to  the  blade 
may  occur  predominantly  in  the  radial  direction.  As  was  done  previously, 
the  rotor  is  arsumed  to  continue  to  translate  in  level,  constant-speed 
flight  during  the  short  time  interval  for  which  transient  responses  are 
computed.  Since  the  method  requires  that  the  wake  geometry  be  specified, 
the  assumption  is  made  that  the  vertical  displacement  of  the  wake  results 
from  the  induced  velocity  yielded  by  momentum  theory,  as  determined  by 
the  instantaneous  thrust  of  the  rotor.  The  aerodynamic  forces  and 
moments  are  based  on  quasi-static  considerations  and,  therefore,  the 
omission  of  shed  vorticity  in  the  wake  is  believed  to  be  consistent  with 
the  adopted  quasi -static  approach. 

Loadings,  bending  moments,  and  modal  responses  are  obtained 
for  a  teetering  rotor  experiencing  both  a  100-  ana  a  50 -percent  power  loss. 
Similarly,  responses  are  computed  for  a  "rigid"  rotor  experiencing  a 
100-percent  power  lobs  and  for  a  "rijid"  rotor  being  braked  from  an  initial 
operational  condition  corresponding  to  100-percent  rpm  and  a  nominal 
value  of  zero  lift.  This  "rigid"  rotor  is  also  started  up  from  a  stopped 
position.  Examination  of  the  results  obtained  in  these  computations 
indicates  (1)  that  reasonable  aerodynamic  loadings  can  be  predicted 
for  all  values  of  the  radial  to  tangential  velocity  ratio  and  (2)  that  there 
exist  no  loading  or  structural  response  phenomena  for  the  examined 
flight  condition  which  would  indicate  that  power  loss,  braking,  and 
rotor  stopping  have  serious  consequences  with  respect  to  the  creation 
of  unusual  blade -wake  interaction  phenomena. 
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INTRODUCTION 


Accurate  prediction  of  blade  loadings,  motion,  and  stress  experi¬ 
enced  by  helicopter  rotors  in  forward  flight  requires  that  a  multitude  of 
aerodynamic  and  structural  dynamic  phenomena  be  modeled  in  a  compre¬ 
hensive  manner.  This  analysis  objective  has  led  to  the  performance  of 
many  studies  [References  1  through  8]  which  have  demonstrated  that 
mathematical  models  can  provide  considerable  assistance  in  isolating 
the  factors  responsible  for  the  vibratory  behavior  exhibited  by  rotary 
wings.  In  particular,  the  development  of  means  for  calculating  the 
temporal  and  spacial  distribution  of  inflow  induced  by  the  rotor  wake 
[References  9  through  111  must  be  cited  as  a  major  forward  step.  This 
development  has  made  it  possible  to  eliminate  the  classical  simplifying 
assumption  of  uniform  inflow  from  rotor -loading  calculations,  which  is  a 
desirable  procedure  when  there  is  need  to  predict  the  higher  harmonic 
components  of  the  excitation  and  response  of  the  dynamic  system  repre¬ 
sented  by  the  rotary  wing. 

Wake  models  enabling  the  calculation  of  nonuniform  inflow  were 
applied  initially  in  predicting  the  cyclic  or  periodic  behavior  of  rotary 
wings  in  steady,  forward  flight.  More  recently,  under  the  sponsorship 
of  USAAVLABS,  a  study  [Reference  12]  was  conducted  for  the  purpose 
of  developing  a  computer  methodology  suitable  for  obtaining  the  non¬ 
periodic  response  of  a  rotary  wing  caused  by  a  transient  input  of  blade 
pitch.  In  this  latter  study,  it  was  assumed  that  the  forward  flight  tra¬ 
jectory  of  the  rotor  hub  does  not  vary  during  the  time  interval  of  interest. 
Further,  blade  flexibility  was  neglected  since  the  primary  objective  was 
to  develop  a  practical  definition  of  the  transient  vorticity  and  geometry 
of  a  rotor  wake  to  demonstrate  the  feasibility  of  computing  the  air  loads 
produced  during  the  transient  interval.  Accordingly,  a  mathematical 
model  was  generated  to  describe  a  rotor  rotating  at  a  fixed  angular 
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velocity  and  possessing  only  a  flapping  degree  of  freedom.  As  was 
done  in  studies  of  steady* state  blade  loads,  the  transport  velocities 
of  elements  of  wake  vorticity  were  estimated  and  specified,  a  priori. 
A  subsequent  comparison  of  calculations  and  experimental  data  dem¬ 
onstrated  that  the  prediction  of  nonperiodic  air  loadings  on  rotary 
wings  is  a  feasible  task. 

The  successful  completion  of  this  preliminary  study  led  to  a 
recommendation  for  introducing  additional  degrees  of  freedom  into 
the  model  such  that  blade  loadings  and  elastic  motions  in  the  flapwise, 
edgewise, and  pitchwise  directions,  as  well  as  the  angular -velocity 
of  the  rotor  shaft,  could  be  predicted  as  functions  of  transient  inputs 
of  blade  pitch  and/or  shaft  torque.  The  recommended  research 
objective  was  to  obtain  an  understanding  of  the  dynamics  of  rotor 
blades  following  a  power  failure,  for  example,  and  to  investigate 
further  rotor  behavior  during  starting  and  stopping  conditions  such 
as  would  occur  in  stopped-rotor  designs.  Specifically,  it  was  recom¬ 
mended  that  a  study  be  made  to  determine  the  transient  response 
characteristics  of  rotary  wings  as  a  function  of  (1)  the  blade-wake 
relationships  and  (2)  the  inertial,  elastic,  and  aerodynamic  coupling 
forces  and  moments  that  prevail  during  the  transient  interval. 

In  making  this  recommendation,  it  was  recognized  that  the 
high  advance  ratios,  associated  with  stopped-rotor  designs,  cause 
flow  conditions  in  which  the  existing  aerodynamic  analysis  procedures 
for  rotors  are  not  valid.  The  development  and  validation  of  theories 
suitable  for  predicting  aerodynamic  loads  on  rotor  blades  at  high 
advance  ratios  are  held  to  be  a  major  research  task,  whose  execution 
has  recently  been  initiated.  Nevertheless,  it  was  essential  that  some 
effort  be  made,  as  part  of  the  study  documented  in  this  report,  to 
come  to  grips  with  this  problem.  This  was  done.  From  the  research 
point  of  view,  it  is  believed  that  the  aerodynamic  model  developed  in 
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this  study  constitutes  a  significant  first  step  toward  the  resolution  of 
the  analysis  difficulties  posed  by  the  operation  of  rotors  at  very  high 
advance  ratios. 

Accordingly,  the  first  part  of  this  report  is  devoted  to  pre¬ 
senting  the  reasoning  behind  the  procedure  adopted  to  compute  loadings 
in  the  presence  of  large  components  of  radial  flow.  Information  of 
interest  to  the  general  reader  ie  given  in  the  body  of  the  report,  and 
a  more  detailed  accounting  of  certain  facets  of  the  aerodynamic  mocel 
are  given  in  Appendix  II.  Subsequent  to  the  discussion  of  the  develop¬ 
ment  of  an  aerodynamic  theory  having  increased  validity  at  high  advanc 
ratios,  a  brief  description  is  given  of  the  overall  simulation  model, 
again  with  the  details  being  placed  in  appropriate  appendixes.  The 
computational  methodology  employed  to  obtain  numerical  results  is 
also  described  briefly.  This  section  is  followed  b/  a  discussion  of  the 
results  of  calculations  that  were  made  using  the  developed  model.  Con 
elusions  drawn  from  these  results  and  recommendations  for  further 
research  conclude  the  report. 

Throughout  the  report,  TBL-I  refers  to  the  computational 
model  developed  and  reported  in  Reference  [12]  for  Transient  Blade 
Loads,  with  TBL-II  used  to  designate  the  computational  model  devel¬ 
oped  in  the  current  vtudy. 
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STATEMENT  OF  THE  PROBLEM 


The  examination  of  blade  loadings,  motions,  and  stresses 
occurring  during  a  power  loss  or  during  a  gradual  starting  or  stopping 
process  is  facilitated  by  expanding  the  analysis  developed  in  Refer¬ 
ence  [12]  to  determine 

(1)  the  inplane  loads 

(2)  the  pitching  momenta  about  the  feathering  axis 
(i.  e. ,  the  control  loads) 

(3)  the  rigid -blade  motions  (flap,  lead -lag,  pitch) 

(4)  the  flexible  blade  motions  (flapwise  and  chordwise 
bending,  blade  twist)  and 

(5)  the  angular -velocity  response  of  the  rotor  shaft 

as  caused  by  time -varying  changes  in  swashplate  displacement  and 
shaft  torque. 

If  a  rotor  is  subjected  to  a  negative  angular  acceleration  as 
a  result  of  a  power  loss  or  a  braking  torque,  advance  ratios  of  increasing 
magnitude  are  encountered,  which,  in  the  limit,  reach  infinity  when  the 
rotor  stops.  Thus,  the  rotor  blades  will  be  exposed  to  a  variety  of 
flow  geometries  that  cannot  be  described  adequately  by  existing  aero¬ 
dynamic  analyses.  It  follows  that  the  requirement  tc  calculate  aero¬ 
dynamic  loadings  at  high  advance  ratios  constitutes  the  major  problem 
to  be  resolved  in  this  investigation.  In  prior  studies  of  rotor  air  loads, 
the  radial  components  of  flow  (i.  e.  ,  flow  along  the  blade)  have  been 
neglected,  thereby  permitting  the  use  of  two-dimensional  section  data 
for  lift,  drag,  and  pitching  moment.  Discrepancies  between  calculated 
and  measured  performance  at  advance  ratios  in  excess  of  0.  3  (e.  g. , 
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References  [13]  and  [14])  have  indicated  the  inadequacy  of  two- 
dimensional  theory  at  high  advance  ratios.  The  question  arises  as 
to  the  upper  bound  of  advance  ratio  for  which  a  two-dimensional 
treatment  is  satisfactory. 

In  addition  to  the  computational  complexities  introduced  by 
the  use  of  a  more  elaborate  aerodynamic  analysis,  the  inclusion  of  the 
lead-lag,  feathering,  and  blade  flexibility  degrees  of  freedom  con¬ 
stitutes  a  significant  expansion  of  the  computing  program  used  in 
TBL-I.  A  further  complication  encountered  in  TBL-II  is  the  loss  of 
the  linear  relationship  between  blade  azimuth  and  time  as  a  result  of 
shaft  torque  being  treated  as  one  of  the  independent  variables  of  the 
study. 


In  retrospect,  it  may  be  stated  that  the  extension  of  the 
analysis  procedure  described  in  Reference  [12]  has  consisted  of 
three  major  tasks.  The  first  task  was  the  development  of  equations 
of  motion  applicable  to  a  variable-rpm  rotor  (i.  e. ,  the  rotor  has 
shaft  displacement  as  a  degree  of  freedom)  in  which  each  blade  can 
be  assumed  to  have  as  many  as  11  degrees  of  freedom  (three  rigid- 
blade  displacements,  six  coupled  flapwise  and  chordwise  bending 
modes,  and  two  uncoupled  torsion  modes).  The  second  task  was  the 
generation  of  an  aerodynamic  model  (namely,  a  system  of  horseshoe 
vortices  and  control  points)  that  yields  a  distribution  of  bound  vorticity 
which  can  be  converted  to  representative  loadings  for  any  ratio  of 
radial  to  tangential  flow  velocities.  The  third  task  consisted  of  gen¬ 
erating  a  machine -computing  program  corresponding  to  the  developed 
mathematical  model. 
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THE  DEVELOPMENT  OF  AN  AERODYNAMIC  MODEL  SUITABLE 


FOR  CALCULATING  LOADINGS  AT  HIGH  ADVANCE  RATIOS 


The  aerodynamic  model  used  in  TBL-I  contained  the  following 
assumptions: 

1.  The  wake  is  trailed  from  the  trailing  edge  of  the 
blade  and  is  represented  by  discrete  vortex  lines 
trailed  as  a  consequence  of  a  stepwise  rate  of 
change  of  the  total  circulation  along  the  span. 

2.  Each  segment  of  the  blade  is  treated  in  accord¬ 
ance  with  the  Prandtl  assumption  that  the  bound 
vorticity  on  a  given  section  extends  continuously 
in  the  spanwise  direction  resulting  in  a  chordwise 
distribution  of  bound  vorticity  that  is  characteristic 
of  the  two-dimensional  airfoil.  The  consequence  of 
this  simplification  is  that  replacement  of  the  chord- 
wise  distribution  of  bound  vorticity  by  a  lifting  line 
at  the  quarter  chord  results  in  the  requirement  for 
flow  tangency  at  the  three-quarter  chord  being 
equivalent  to  finding  the  lift  on  a  two-dimensional 
section, operating  at  an  effective  angle  of  attack 
equal  to  the  free-stream  angle  of  attack  minus  the 
angle  of  attack  induced  by  the  free  wake. 

By  restating  the  above  assumptions  in  mathematical  form,  the 
effective  "strip"  nature  of  the  analysis  used  in  TBL-1  becomes  clear. 
With  the  chordwise  distribution  of  vorticity  corresponding  to  a  two- 
dimensional  flow  replaced  by  a  "lumped"  bound  vortex  at  the  quarter 
chord  (extending  to  plus  and  minus  infinity),  the  velocity  induced  by 
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this  bound  vortex  at  the  three-quarter  chord  is 


n 

2rrbj 


(1) 


where 

uf}  is  the  velocity  induced  by  the  bound  vortex  of 
infinite  length 

/ j  is  the  strength  of  the  bound  vortex  at  the  jth 

span  location  with  this  strength  assumed  to 
extend  to  plus  and  minus  infinity 

bj  is  the  half  chord  at  the  j **  span  location. 

For  the  assumed  two-dimensional  flow,  the  requirement  for 
flow  tangency  at  the  three-quarter  chord  may  be  stated  as 

<"b  +  (j)  +  (af  Uj)^  -  0  (2) 


where 

is  the  velocity  induced  by  the  free  vorticity  in  the  wake 
at  the  three-quarter  chord  of  the  j **  span  segment 

and 

**  comPonent  free-stream  velocity  per¬ 

pendicular  to  the  blade  at  the  three-quarter  chord  where  a9 
and  U  are  defined  only  in  terms  of  velocities  lying  in  a 
plane  perpendicular  to  the  span  of  the  blade. 

Equations  (1)  and  (2)  yield  that 

Seiilj.  (3) 
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On  identifying  the  constant,  Zrr ,  as  the  theoretical  lift-curve 
slope  for  a  two-dimensional  flat  plate,  it  is  seen  that  the  following  holds 
for  a  real  airfoil: 

ri  ■  h  Uj  CtU) 

where 


Cj(j)  is  an  experimentally  determined  lift  coefficient  for  a 
two-dimensional  section  at  an  effective  angle  of  attack  equal 
to 

'  <4- 

It  is  seen  that  the  above  model  contains  the  full  influence  of 
the  three-dimensional  wake  and,  furthermore,  retains  the  advantage 
of  using  two-dimensional  section  characteristics  for  the  airfoil.  Although 
it  is  possible  to  compute  an  asymmetrical  loading  increment  that  results 
from  a  radial  component  of  the  free-stream  velocity  reacting  with  chord- 
wise  components  of  the  bound  vorticity,  this  computation  is  approximate, 
at  best,  since  the  boundary  condition  as  applied  does  not  include  this  radial 
component  of  the  free  stream. 


In  order  to  obtain  a  calculation  procedure  that  is  applicable  to 
the  flow  regime  wherein  a  high -aspect-ratio  wing  becomes  a  very  low- 
aspect-ratio  lifting  surface,  the  foregoing  "strip"  analysis  must  be 
modified.  One  obvious  alternative  is  to  utilize  a  lifting-surface  approach. 
Such  an  approach  would  involve  a  large  number  of  control  points  and  would 
necessitate  computing  the  velacity  induced  both  by  the  bound  vorticity 
on  the  blade  and  the  free  vorticity  in  the  wake  at  each  and  every  control 
point.  For  very  practical  reasons,  namely,  computer  capacity,  it  was 
necessary  that  the  number  of  control  points  be  restricted  to  a  minimum. 

In  addition,  it  was  necessary  to  adopt  a  system  of  vortices  for  repre¬ 
senting  a  wing  in  the  presence  of  oblique  flow  such  that  the  new  system 
reduces  to  a  vortex  model  that  is  appropriate  fcr  computing  loadings  at 
low  advance  ratios  and  is,  in  practice,  equivalent  to  that  used  previously 
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for  the  "strip"  analysis.  A  means  for  checking  the  reasonableness 
of  the  loadings  yielded  at  low  advance  ratios  by  the  modified  aero¬ 
dynamic  calculation  procedure  is  thereby  assured. 

During  the  process  of  developing  an  aerodynamic  model  suitable 
for  TBL-II,  an  examination  of  the  literature  indicated  that  moderately 
low-aspect-ratio  lifting  surfaces  are  characterized  by  a  relatively 
constant  distribution  of  loading  over  the  span  (References  [15]  through 
[17]).  Moreover,  this  relatively  constant  spanwise  loading  tends  to 
prevail  at  all  positions  along  the  chord.  Further,  a  sheet  of  vorticity 
emanates  from  the  side  edges  of  the  low-aspect-ratio  surface  and  this 
sheet  proceeds  to  roll  up  into  two  pronounced  wing-tip  vortices  lying 
above  the  plane  of  the  wing.  For  a  very  low -aspect-ratio  surface,  it 
appears  appropriate  to  approximate  the  observed  constant  loadings 
across  the  span  and  the  existence  of  the  side -edge  vortex  sheet  by  the 
system  of  horseshoe  vortices  diagrammed  in  Figure  1. 


Figure  1.  HORSESHOE-VORTEX  REPRESENTATION  OF  A  10W-ASPECT-RATI0 

SURFACE  \.ITH  SHARP  EDGES  AT  THE  WINS  TIPS  AT  ANGLE  OF  ATTACK,*. 
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Figure  1,  me  drawn,  is  representative  of  a  rotor  blade  exposed 
to  a  completely  radial  flow.  The  sharply  defined  side -edge  separation 
would  exist  for  a  lifting  surface  having  very  sharp  edges  at  the  wing 
tips.  On  a  helicopter  rotor  blade,  these  side  edges  consist,  however, 
of  a  rounded  leading  edge  and  a  sharp  trailing  edge  where  "leading" 
and  "trailing"  are  defined  in  terms  of  a  flow  that  is  primarily  tangen¬ 
tial  to  the  blade.  With  these  edge  properties  present,  questions  arise 
as  to  the  manner  in  which  side -edge  separation  would  develop  on  a 
rotor  blade  exposed  to  an  exclusively  radial  flow.  This  problem  of 
defining  the  geometry  and  location  of  the  vortex  sheets  emanating  from 
a  lifting  surface  becomes  even  more  difficult  when  a  typical  rotor  blade 
is  exposed  to  a  flow  which  is  other  than  completely  radial  or  is  exposed 
to  a  flow  of  varying  obliquity  as  caused  by  a  slow  rotation  rate. 

In  view  of  the  large  gap  in  our  understanding  of  the  manner 
in  which  a  rounded  edge  delays  the  creation  of  the  vortex  sheet  asso¬ 
ciated  with  the  long  side -edge  dimension,  the  assumption  is  made,  in 
the  TBL-II  analysis,  that  vorticity  is  never  released  into  the  free  stream 
simultaneously  from  the  two  radial  edges  of  the  blade.  Instead,  it  is 
assumed  that  vorticity  is  released  at  the  sharp  trailing  edge  of  the  blade 
whenever  the  local  component  of  tangential  flow  is  positive  and  from 
the  rounded  leading  edge  whenever  the  local  component  of  tangential 
flow  is  negative.  In  this  manner,  the  wake  geometry  existing  at  low 
advance  ratios  is  preserved,  such  that  the  single  vortex  sheet  asso¬ 
ciated  with  the  radial  distribution  of  loading  in  predominantly  tangential 
flow  becomes  a  single  vortex  sheet  emanating  from  a  long  side  edge 
when  the  flow  is  exclusively  radial.  Thus,  the  wing-tip  separation, 
which  is  produced  by  two  sharp  side  edges  of  a  low-aspect-ratio  sur¬ 
face  exposed  to  a  radial  flow  (see  Figure  1),  will  take  the  form  diagrammed 
in  Figure  2(a)  as  a  result  of  the  "single -sheet  assumption".  Further, 
for  an  oblique  flow  condition  ( e.  g. ,  for  the  case  of  the  blade  located  in 
the  retreating  half  of  the  rotor),  it  is  assumed  that  the  bound  and  adjacent 
free  vorticity  would  have  the  geometry  diagrammed  in  Figure  2(b). 
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(•) 


A  ROTOR  BLADE  EXPOSEO  TO 
A  RADIAL  FLOW;  u  =  o°. 


b  A  RETREATING  ROTOR  BLADE  WITH 
THE  INBOARD  SECTION  LOCATED  IN 
A  REVERSED-FLOW  REGION. 


Figure  2.  ASSUMED  HORSESHOE- VORTEX  REPRESENTATION 


In  this  manner,  it  is  possible  to  extend  the  "strip"  analysis  to 
include  the  regime  of  highly  oblique  flow  in  which  a  high-aspect-ratio 
wing  becomes  a  very  low-aspect-ratio  lifting  surface.  First,  it  is 
assumed  that  the  high-aspect-ratio  loadings,  as  yielded  by  the  model 
diagrammed  in  Figure  3,  can  be  obtained  with  little  loss  in  accuracy 
by  shifting  the  spanwise  bound  vortices  from  the  quarter  chord  to  the 
leading  edge  and  placing  the  control  point  at  the  midchord  (see  Figure  4). 


Figure  3.  HORSESHOE-VORTEX  REPRESENTATION  OF  A  ROTARY  Ml  NO 
AT  LOW  ADVANCE  RATIOS. 
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The  geometry  diagrammed  in  Figure  4  (specifically, with  the  bound 
element  of  spanwise  vorticity  located  at  the  leading  edge)  prevails 
whenever  a  given  blade  section  is  in  a  normal  flow  region.  The 
geometry  diagrammed  in  Figure  4  is  reversed  (i.  e. ,  the  bound  ele¬ 
ment  of  spanwise-oriented  vorticity  is  located  at  the  trailing  edge) 
whenever  a  given  blade  section  is  located  in  the  reversed  flow  region 
(i.  e. ,  the  local  tangential  component  of  flow  is  negative). 


Figure  4.  MODIFIED  HORSESHOE- VORTEX  REPRESENTATION  OF  A  ROTARY  WING: 
BOUND  VORTEX  AT  LEADING  EDGE;  CONTROL  POINT  AT  MIDCHORD. 


In  addition  to  the  assumed  shift  in  (1)  the  location  of  the 
radially  oriented  bound  vortex  and  (2)  the  point  at  which  the  condition 
of  flow  tangency  is  to  be  satisfied,  it  is  apparent  that  the  so-called 
Prandtl  simplification  cannot  be  used  to  satisfy  the  boundary  condition 
on  the  low-aspect-ratio  blade.  Elimination  of  the  Prandtl  simplification 
means  that  it  is  necessary  to  compute,  at  each  control  point,  the  veloc¬ 
ities  induced  by  each  radial-  and  chordwise-oriented  element  of  vor¬ 
ticity  assumed  to  represent  the  lifting  surface  as  well  as  the  velocities 


12 


induced  by  the  free  wake.  Further,  the  boundary  condition  must  be 
generalized  to  include  all  components  of  the  total  flow  vector  relative 
to  the  blade  segment  under  consideration.  The  equations  that  yield 
the  unknown  circulation  strengths  on  the  blade  accordingly  derive 
from  a  "three-dimensional"  boundary  condition;  viz. , 

*  “^(j)  +  Oc(j)  *  0  (4) 

where 

U^(j)  is  the  component  of  the  total  free-stream  /elocity 
perpendicular  to  the  blade  at  the  control  point  of  the 
radial  segment 

and 

or^(j)  and  u^r(J)  can  be  expressed  symbolically  as 
( j)  » 

&~<j)  •  Z.  EL  (  ) 

where 

5= v  and  <7^m^are  the  induced  velocity  coefficients  of 
the  bound  and  free  vorticity,  respectively. 

Prior  to  accepting  the  above -outlined  aerodynamic  model 
as  a  means  for  determining  loadings  at  high  advance  ratios  and  at 
azimuth  angles  in  the  vicinity  of  0  and  180  degrees,  computations 
were  performed  to  ascertain  whether  Equations  (4)  and  (5)  yield 
loadings  in  agreement  with  low-aspect-ratio  wing  theory;  namely, 
the  linear  theory  applicable  to  vanishingly  small  angles  of  attack. 

If  the  angle  of  attack  of  the  blade  diagrammed  in  Figure  2(a)  is 
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assumed  to  approach  sero,  the  t.  ie  vortices  trailing  radially  in  the 
free  airstream  will  lie  in  the  plane  of  the  wing  and  the  resulting 
system  of  vortices  reduces  to  the  system  diagrammed  in  Figure  5. 


Figure  5.  HORSESHOE- VORTEX  REPRESENTATION  OF  LOW-ASPECT-RATIO  WINS 
AT  AN  INFINITESIMAL  ANGLE  OF  ATTACK. 


A  loading  calculation,  equivalent  to  linear  airfoil  theory,  consists  of 
evaluating  the  induced  velocity  coefficient  matrix  corresponding  to  the 
system  of  vortices  diagrammed  in  Figure  5.  An  aspect  ratio  of  0. 1 
was  selected.  Calculations  were  performed  for  both  10  and  20  radial 
segments  involving  10  and  20  control  points,  respectively. 
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For  the  assumed  radial  flow  with  no  sideslip  component  present, 
the  boundary  condition  expressed  for  all  control  points  has  the  form 


l-b 

The  solution  of  the  above  set  of  equations  for  the  unknown  circulations,  /J  , 
can  be  written  symbolically  as 


’ 

- 

-/ 

’  -  / 

n 

/T  ■ 

M 

c/ 

.  •  4 

( ocU\ 

: 


The  quantities,  /J /(Jot*  are  "lumped"  values  of  the  chordwise  distribution 
of  loading  and  can  be  converted  to  a  continuous  loading  distribution  by 
obtaining  the  nondimensional  variable,  where  x  is  the  distance 

along  the  chord  and  f  is  the  local  strength  of  the  spanwise  vorticity  as 
distributed  along  the  chord.  The  results  of  the  calculation  are  presented 
in  Figure  6, and  it  is  observed  that  the  chordwise  distribution  of  circu¬ 
lation  is  characteristic  of  that  exhibited  by  very  low-aspect-ratio  surfaces. 


Figure  6.  CHORDWISE  DISTRIBUTION  OF  CIRCULATION  ON  A 
RECTANGULAR  WING  OF  0.10  ASPECT  RATIO. 
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On  noting  that  the  total  lift  is  given  by  the  following  identity: 

[c^lL)dx  .  su'Cibs 

where 

b  is  the  half  chord 
j  is  the  span 

and 

■y  is  the  local  circulation  with  constant  strength  across 
the  span, 

we  find  that  Cl/<*  is 

a,  b  J0  ax 

r 

On  performing  the  integration,  we  obtain  z^.  =  0.  209.  If  the  spanwise 
distribution  of  circulation  had  been  assumed  to  be  elliptical  instead  of 
rectangular,  the  lift-curve  slope  would  be  0.  209  x  1L  =  .  Ib4,  which 
compares  very  favorably  with  the  value  of  0.  157  yielded  by  the  low- 
aspect-ratio  wing  theory  of  Jones  (Reference  [  18]). 

The  above-described  system  of  horseshoe  vortices  and  locations 
for  satisfying  the  condition  of  flow  tangency  is  thus  found  to  yield  accept¬ 
able  chordwise  loadings  for  the  low-aspect-ratio  case  and  is  known  to 
be  a  satisfactory  approximation  of  the  spanwise  distribution  of  loading 
for  high-aspect-ratio  surfaces.  As  a  result  of  the  assumptions  made 
with  respect  to  the  convection  of  the  vorticity  trailed  into  the  air  stream, 
this  aerodynamic  model  contains  the  mechanism  wherein  the  loading  of 
the  low-aspect-ratio  surface  becomes  a  nonlinear  function  of  the  angle 
of  attack.  Demonstration  of  the  adequacy  of  the  model  at  the  two  extremes 
of  aspect  ratio  forms  the  oasis  of  the  heuristic  argument  that  the  model 
will  yield  adequate  results  throughout  the  spectrum  of  effective  aspect 
ratios  and  advance  ratios  encountered  by  the  blades  of  a  rotary  wing. 
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Equation  (5)  shows  that  the  simple  relationship  between  u^(j) 
and  /j  given  by  Equation  (1)  does  not  hold  for  the  three-dimensional 
definition  of  the  boundary  condition.  Accordingly,  it  is  not  possible 
to  apply  the  three-dimensional  boundary  condition  and  obtain  a  result 
similar  to  Equation  (3)  which  can  be  interpreted  in  terms  of  empirical 
data  obtained  from  two-dimensional  tests  of  various  airfoil  sections. 
Note  also  that  the  horseshoe-vortex  model  in  the  proposed  form  will 
yield  loadings  that  derive  only  from  circulatory  forces.  However,  it 
is  most  desirable  that  the  model  include  means  for  approximating 
the  influence  of  viscosity  on  the  slope  of  the  lift  curve  and  means 
for  introducing  forces  resulting  from  viscous  effects.  Further,  the 
influence  of  Mach  number  on  the  two-dimensional  lift-curve  slope 
should  be  accounted  for  in  this  analysis  in  order  that  the  model  produce 
the  same  results  at  low  advance  ratios  as  would  be  obtained  from 
TBL-1.  Finally,  it  was  desirable  that  means  be  devised  for  predicting 
the  center  of  pressure  location  along  the  tangential  direction,  irre¬ 
spective  of  whether  the  flow  past  the  blade  was  primarily  tangential 
or  radial,  in  order  to  compute  pitching  moments  about  the  feathering 
and  elastic  axes.  These  objectives  were  achieved  by  appropriate 
extensions  and  modifications  of  the  vortex  model,  the  details  of  which 
are  described  in  Appendix  I.  The  logic  for  these  modifications  was 
dictated  by  the  requirement  that  the  application  of  the  proper  boundary 
conditions, when  the  flow  is  primarily  two-dimensional,  should  reduce 
the  model  to  the  one  appropriate  for  the  high-aspect-ratio  analysis 
as  used  in  TBL-I.  Appendix  I  also  describes  the  manner  in  which 
the  system  of  horseshoe  vortices  placed  at  the  blade  are  altered 
whenever  a  portion  of  the  blade  is  in  the  reversed  flow  region. 

In  developing  the  load-calculation  procedure  outlined  above, 
the  question  arose  as  to  the  manner  in  which  to  introduce  the  influence 
of  a  local  stalled  flow  into  the  three -dimensional  boundary  condition. 
Examination  of  Equation  (4)  shows  that  this  formulation  of  the  boundary 
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condition  has  no  built-in  "mechanism"  for  limiting  the  circulation 
on  the  blade  as  might  be  caused  by  flow  separation.  Since  it  appeared 
desirable  that  some  form  of  stall  limit  be  included,  at  least  for  the 
case  where  the  flow  is  primarily  two-dimensional,  means  are  provided 
for  limiting  the  bound  circulation  computed  at  low  advance  ratios  to 
the  value  given  by  the  maximum  lift  coefficient  obtained  in  two-dimen¬ 
sional  wind-tunnel  tests. 


This  limiting  procedure  is  readily  implemented  by  defining 
the  maximum  circulation  in  the  j**1  horseshoe  vortex  to  be  given  by 


r ) 

J/m»u 


(6) 


where 


and 


U L  ,  Uf  ,  and  Ut  are  the  tangential,  normal,  and  radial 
"i  V  V  x/ 

components  of  the  free-stream  velocity  at  the  j*"  radial 

segment 


Cj(/mAXia  the  local  maximum  lift  coefficient  produced 
in  two-dimensional  flow. 


Note  that  at  low  values  of  advance  ratio,  Equation  (6)  yields 


~  kj  Uff. 


(7) 


and  for  the  flow  condition  with  predominantly  radial  flow,  yields 


rj)ma.X  ~  • 


(8) 


It  is  necessary  to  show  that,  whereas  Equation  (7)  does  im¬ 
pose  a  limit  on  Cj  for  a  two-dimensional  flow,  the  same  two-dimensional 
value  of  Cj(j)  does  not  impose  a  similar  limit  when  the  flow  is  predom¬ 
inantly  radial.  Fr  r  the  radial  flow  case,  the  peak  value  of  circulation 


. 


* 

i 

occurs  at  the  leading  edge  (i.  e.  ,  at  the  tip  of  the  blade)  and  for  the 
0.  10  aspect-ratio  surface  treated  earlier  is  given  by 

-  '■*“,<*  • 

On  defining 

“  u(t  &)snu 

and  on  using  Equation  (8),  the  stall  angle  of  attack,  under  conditions 
of  radial  flow,  is 

On  assuming  -  1.  3  (a  reasonable  assumption  for  low  Mach 

numbers),  the  stall  angle  of  attack  for  the  radial-flow  case  is  approxi¬ 
mately  37  degrees.  This  angle  is  sufficiently  large  to  make  it  apparent 
that  a  limitation, based  on  the  maximum  lift -coefficients  produced  by 
the  airfoil  sections  in  two-dimensional  flow,  does  not  impose  a  sig¬ 
nificant  restriction  on  the  values  of  lumped  circulation  that  define 
the  loading  when  the  blade  is  exposed  to  a  predominantly  radial  flow. 
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THE  DYNAMICS  OF  A  ROTOR-WAKE  SYSTEM 


A  discussion  of  the  simulation  developed  to  represent  the  time- 
varying  response  of  a  flexible,  rotary  wing  is  facilitated  by  considering 
the  various  simulation  submodels  that  combine  to  represent  the  total 
dynamic  system.  These  submodels  <ne  arbitrarily  designated  here 
as:  (1)  the  "wake  model",  (2)  the  "aerodynamic  model",  and  (3)  the 
"structural  model".  Their  general  character  and  interaction  with 
each  other  are  discussed  in  this  section  together  with  the  assumptions 
and  simplifications  that  were  made  for  purposes  of  achieving  a  tractable 
computer  program.  A  more  complete  description  of  the  submodels  will 
be  found  in  appendixes  to  this  report. 

It  proved  convenient  to  develop  the  TBL-II  simulation  with  two 
separate  aerodynamic  submodels.  One  of  the  aerodynamic  models  is  a 
slight  extension  of  that  employed  in  TBL-I  in  that  two-dimensional  airfoil 
data  are  used  to  obtain  both  lift  and  drag  loadings  on  the  blade.  The 
second  model  is  that  devised  to  calculate  the  loadings  on  blades  exposed 
to  flow  of  high  obliquity.  The  motivation  for  incorporating  an  aerodynamic 
"strip"  theory  into  TBL-II  stemmed  from  the  need  to  obtain  a  convenient 
checking  procedure  for  the  digital  program  itself  and  from  the  need  to 
ascertain  whether  the  aerodynamic  model  suitable  for  high-advance-ratio 
calculations  produces  correct  results  at  low  advance  ratios. 

In  developing  the  TBL-II  simulation,  a  continual  compromise 
had  to  be  made  between  the  objective  of  achieving  broad  generality  and 
the  requirement  for  remaining  within  the  capacity  of  the  available  com¬ 
puting  equipment.  In  practice,  it  is  very  difficult  to  make  this  choice 
since  it  is  not  always  clear  whether  a  particular  feature  of  the  model 
is  a  significant  determinant  of  the  behavior  to  be  exhibited  by  the  overall 
system.  Experience  and  judgment  are  called  upon  in  the  selection  of 
the  model  features  outlined  in  the  following  sections. 
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AXES,  ROTOR  GEOMETRY,  KINEMATICS 

Figure  7  shows  the  x-y-f  axis  system  selected  to  define  the 
geometry  of  the  wake.  Note  that  the  origin  is  placed  at  the  hub  of  the 
rotor,  with  the  x-axis  directed  rearward,  parallel  to  the  horizontal 
direction  of  flight.  Figure  7  shows  that  the  rotor  shaft  is  assumed  to 
tilt  only  in  the  x-j  plane.  Thus,  we  ignore  any  lateral  tilt  of  the  rotor 
shaft,  since  the  tilt  angle  is  usually  small  in  trimmed,  level  flight.  Both 
a  nonrotating  and  a  rotating  set  oi  axes  (with  the  axis  aligned  along  the 
rotor  shaft)  are  used  to  define  the  azimuthal  displacement  of  the  rotor. 
(See  Figure  7.)  Additional  axes  are  used,  as  shown  in  Figures  8  and  9 


Figure  7.  AXES  USED  TO  DEFINE  WAKE  GEOMETRY  AND  THE 
ROTATION  OF  THE  ROTOR  SHAFT. 
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The  structural  degrees  of  freedom  are  defined  in  terms  of  linear 

r 

displacements  of  the  elastic  axis  in  the  local  chordwise  and  normal  direc¬ 
tions  and  of  angular  displacement  about  the  elastic  axis.  These  displace¬ 
ments  are  designated  in  Figure  9,  for  a  given  radial  location  on  the  blade, 


Figure  9.  DEFINITION  OF  THE  STRUCTURAL  DEGREES  OF  FREEDOM  AND**) 

AND  THE  FEATHERING  DEGREE  OF  FREEDOM  (0). 
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as  i/;<i^and  pt  respectively.  Without  any  deformation  of  the  blade,  the 
elastic  axis  is  assumed  to  be  straight  but  not  necessarily  parallel  to  the 
feathering  axis.  The  location  of  the  elastic  axis  with  respect  to  the 
feathering  axis,  Xy  ,  of  the  blade  is  given  (as  shown)  by  the  dimension  eo  , 
this  dimension  being  a  function  of  radial  position  along  the  blade.  At 
any  given  radial  position,  a  axis  system  is  used  to  define  the  radial, 

chordwise,  and  normal  directions  with  respect  to  a  plane  tangent  to  the 
surface  of  the  blade  at  a  specified  control  point.  Figure  10  shows  the 
definition  used  for  angle  of  attack  in  the  two-dimensional  "strip"  analysis. 


(Jg  ~  AIR  VELOCITY  RELATIVE  TO  THE  BLADE 
f  ALONG  THE  POSITIVE  £  AXIS 


U.  -  AIR  VELOCITY  RELATIVE  TO  THE  BLADE 
7  ALONG  THE  NEGATIVE  n  AXIS 

-/  *4 

ANGLE  OF  ATTACK,  ct  =  to. * 

un 

Figure  10.  DEFINITION  OF  ANGLE  OF  ATTACK  APPLICABLE 
TO  LOW-ADVANCE-RATIO  CALCULATIONS. 


A  slight  inconsistency  exists  in  the  analysis  as  a  result  of  assuming  that 
the  built-in  twist  of  the  blade  is  defined  with  respect  to  the  elastic  axis  in 
deriving  the  structural  equations  of  motion,  whereas  in  the  kinematic 
analysis, it  was  assumed  that  the  built-in  twist  is  defined  with  respect  to 
the  feathering  axis.  In  actual  practice,  a  blade  is  twisted  about  neither 
of  these  axes.  This  discrepancy  is  not  considered  to  be  significant. 
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In  the  three-dimensional  analysis,  is  defined  to  include  the  velocity 
induced  by  the  free  wake  but  does  not  include  the  velocity  components 
caused  by  bound  vorticity.  (Note  that  the  kinematic  analysis  required 
to  express  displacements  and  velocities  in  an  appropriate  reference 
system  is  given  in  Appendix  I. ) 

WAKE  MODEL 

The  wake  model  employed  in  TBL-U  is  similar  to  that  utilized 
in  TBL-I  with  a  few  minor  exceptions.  Although  the  TBL-II  analysis 
contains  higher  frequency  motions  as  a  result  of  simulating  blade 
flexibility,  the  vorticity  shed  parallel  to  the  trailing  edge  of  the  blade, 
as  a  result  of  the  time  rate  of  change  of  bound  vorticity  on  the  blade, 
has  been  neglected  as  was  done  in  TBL-I.  In  effect,  a  quasi-static 
approach  is  used  in  computing  aerodynamic  loadings,  and  this  approach 
is  equivalent  to  stating  that  the  influence  of  the  shed  vorticity  in  the 
wake  is  second  order  in  comparison  with  the  influence  of  the  trailing- 
vortex  system. 

The  indexing  procedure  defined  in  Figure  11  is  identical  to 
that  employed  in  TBL-I  and  serves  as  a  convenient  means  of  identifying 
the  variables  in  the  analysis.  As  before,  m  designates  an  instant  of 
time;  i  ,  a  blade  of  the  rotor;  J  ,  the  spanwise  location  of  a  control 
point;  ib  ,  the  wake  trailed  by  blade  i  ;  m ,  an  azimuthal  position  in 
the  wake;  and  J. ,  a  spanwise  position  in  the  wake  or  the  spanwise 
location  of  a  blade-segment  end  point.  For  computing  the  matrix  of 
induced-velocity  coefficients,  the  Biot-Savart  law  was  applied  as 
described  in  the  appendix  of  Reference  [12]. 

It  was  again  assumed  that  only  the  y  -component  of  induced 
velocity  need  be  retained  in  the  analysis. 
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Figure  11.  INDEXING  PROCEDURE  USED  TO  DESIGNATE  BLADE  SEGMENTS, 
CONTROL  POINTS,  AND  WAKE  LOCATIONS. 


Aa  was  done  in  TBL-I,  the  wake  geometry  is  arbitrarily  spec* 
ified.  However,  the  specification  procedure  has  been  altered.  In 
TBL-H,  the  vertical  transport  velocity  (assumed,  as  before,  to  be  equal 
to  the  momentum  value  of  downwash)  is  continuously  computed  as  a  function 
of  the  instantaneous  value  of  the  total  thrust  produced  by  the  rotor.  This 
modification  has  greatly  simplified  the  computational  procedure  in  that  it 
is  no  longer  necessary  to  predict  in  advance  the  transient  variation  in 
the  mean,  induced  downwash  that  will  occur  as  a  result  of  a  control  in¬ 
put.  In  effect,  the  original  open-loop"  procedure  for  specifying  wake 
geometry  has  been  converted  into  a  "closed-loop"  procedure  based  on 
the  assumption  that  elements  of  trailing  vorticity  are  transported  in  the 
vertical  direction  with  a  constant  velocity  component  equal  to  the  momen¬ 
tum-theory  value  of  the  downwash  that  existed  at  the  instant  that  the  vor¬ 
tex  element  was  trailed  into  the  wake. 
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AERODYNAMIC  MODEL 


For  reasons  of  obtaining  a  numerical  check  at  low  advance 
ratios,  the  computational  program  developed  in  TBL-II  provided  for 
two  aerodynamic  models.  The  first  model  was  that  developed  in 
TBL-I,  subsequently  augmented  to  permit  the  use  of  nonlinear  aero¬ 
dynamic  data  applicable  to  a  360 -degree  range  in  angle  of  attack. 

The  second  model  was  derived  from  the  bound-vortex  and  free-wake 
geometry  that  was  adopted  to  represent  a  lifting  surface  which  alter¬ 
nately  changes  from  a  high-aspect-ratio  to  a  low -aspect-ratio  surface. 

In  both  instances,  it  was  assumed  that  the  continuous  distribution  of 
vorticity  on  the  blade,  with  respect  to  the  radial  coordinate  and  time, 
can  be  adequately  defined  by  using  a  sufficient  number  of  blade  seg¬ 
ments  and  azimuth  intervals,  over  which  segments  and  intervals 
the  vorticity  is  assumed  to  be  constant. 

At  each  solution  interval,  the  unknown  bound  circulations  on 
the  blade  must  have  values  such  that  the  requirement  for  flow  tangency 
is  satisfied  at  each  control  point.  The  equations  for  the  unknown  circu¬ 
lations  accordingly  derive  from  the  boundary  condition.  Since  the 
velocities  induced  at  the  control  points  by  either  the  free  or  bound 
vorticity  can  be  expressed  (by  means  of  the  Biot-Savart  law)  in  terms 
of  the  prevailing  geometry  and  the  strengths  of  all  vortex  elements,  the 
boundary  condition  reduces  to  equations  that  involve  current  and  past 
circulation  strengths  and  the  specified  geometry  of  the  wake. 

In  solving  the  equation  set  resulting  from  the  "strip"  analysis, 
a  nonlinear  variation  of  lift  coefficient  with  angle  of  attack  requires  that 
an  iterative  scheme  be  used.  It  was  convenient  to  start  the  solution 
process  by  first  solving  a  linearized  formulation  of  the  boundary  condition. 

The  linsar  equations  that  result  for  the  "strip"  analysis  (see 
Appendix  in  Reference  [12])  can  be  indicated  symbolically  as 
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where  the  column  matrix,  /  ,  represents  the  unknown  circulation 

strengths.  The  column  matrix,  ^  rjf  •  j  ,  represents  the  portion  of 
the  circulation  at  each  spanwise  segment  determined  by  all  contributions 
to  the  local  angle  of  attack,  other  than  the  contribution  deriving  from  the 
unknown  circulations.  The  square  matrix,  *  *8  derived  from 

induced  velocity  coefficients  that  apply  only  to  elements  of  wake  vorticity 
adjacent  to  the  blade.  These  elements  are  assumed  to  have  the  same 
strength  as  the  unknown  bound  vorticity  on  the  blade. 

The  "geometric"  circulations,  ,  are  determined  by  the 

so-called  "geometric"  angles  of  attack  which  are  functions  of  (1)  the 
orientation  and  velocity  of  a  given  blade  segment  (namely,  blade  kine¬ 
matics)  and  (2)  the  velocity  induced  by  all  elements  of  the  wake,  except 
those  elements  adjacent  to  the  trailing  edge  of  the  blade.  Thus,  the 
downwash  created  by  the  wake  is  arbitrarily  divided  into  two  parts. 

The  first  part  of  the  downwash,  urr  ,  is  given  by  the  following  summation 
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where  the  summation  from  *.i  to  TPi  «  *8LH  represents  the  nonrolled 
wake  exclusive  of  the  elements  adjacent  to  the  blade.  It  should  be  noted 
that  all  of  the  circulations,  r  ,  in  the  above  expression  are  known  since 
they  are  computed  at  an  earlier  time  and  thus  correspond  to  what  can 
be  referred  to  as  the  "old"  wake.  The  second  part  of  the  downwash  is 
designated  as  oSA  and  is  given  by  the  summation 
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where  the  circulations,  f  ,  in  this  expression  are  the  unknown  strengths 
of  the  elements  of  "new"  wake  adjacent  to  the  blade. 

On  obtaining  a  solution  to  the  linearized  equation  set,  wherein 
r9  is  defined  as  a  function  of  ,  the  iterative  process  is  initiated  by 
computing  and  then  correcting  the  angle  of  attack  of  each  blade 
segment  to  reflect  the  influence  of  the  adjacent  ("new")  wake.  A  new 
set  of  lift  coefficients  is  obtained  using  the  experimental  data, and 
the  bound  circulations  of  the  blade  are  redefined.  By  recalculating 
the  downwash  caused  by  the  adjacent  wake,  urA  ,  this  process  is  repeated 
until  the  changes  computed  in  rm-  ■  are  either  less  than  a  specified 
percentage  or  less  than  a  specified  small  number.  (It  should  be  noted 
that  a  Gauss-Seidel  solution  procedure  is  required,  along  with  an 
averaging  technique  for  f  and  c*rA  ,  in  order  to  converge  rapidly  to 
a  stable  solution.  ) 

The  linear  equations  that  result  from  the  "three-dimensional" 
boundary  condition  (see  Appendix  II)  can  be  indicated  symbolically  as 

where  the  column  matrix  j j  j  represents,  as  before,  the 
unknown  circulation  strengths.  The  column  matrix  j  j  repre¬ 

sents  that  portion  of  the  velocity  perpendicular  to  the  blade  (at  the 
control  point)  resulting  from  the  motion  and  displacement  of  the  blade 
and  the  downwash  created  by  the  "old"  wake.  The  square  matrix, 

[cr^'J  ,  represents  the  induced-velocity  coefficients  that  yield  the 
downwash  created  by  all  of  the  unknown  vorticity,  both  bound  and 
free,  where  the  free  vorticity  is  only  that  portion  of  the  wake  imme¬ 
diately  adjacent  to  the  trailing  edge  of  the  blade.  If  particular  values 
of  rm  ■  .  obtained  from  the  solution  of  the  above  equation  set  should 
exceed  a  maximum  value  of  circulation,  calculated  on  the  basis  of  a 
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two-dimensional  stall  condition,  the  vorticity  deposited  in  the  free  wake 
is  arbitrarily  limited  to  the  computed  stall  value.  The  induced  velocity 
components,  SrA  ,  ars  then  computed  using  the  stall-limited  value  of 
circulation. 

The  remaining  features  of  the  aerodynamic  model  are  concerned 
with  the  manner  in  which  the  computed  circulations  for  the  bound  vor¬ 
tices  representing  the  blade  are  converted  to  normal  and  tangential 
loadings  in  pounds  per  foot  of  radial  distance  and  to  pitching  moments 
about  a  radial  axis.  In  the  two-dimensional  "strip"  analysis,  experi¬ 
mental  values  of  lift  and  drag  are  converted  to  a  normal  force  and  a 
chordwise  force.  Further,  the  location  of  the  center  of  pressure  is  ob¬ 
tained  from  available  experimental  data  as  a  combined  function  of  angle 
of  attack  and  Mach  number.  (See  Appendix  II  for  details.  ) 

As  was  discussed  previously,  a  fulfillment  of  the  three-dimen¬ 
sional  boundary  condition  yields  circulatory  forces  that  derive  only 
from  potential  flow  considerations.  Accordingly,  provision  was  made 
for  adding  a  tangential  (i.  e. ,  chordwise,  for  the  high-aspect-ratio 
surface)  viscous  force  to  the  circulatory  forces  which  derive  from 
the  radial  and  tangential  components  of  the  bound  vorticity  representing 
the  blade.  In  accordance  with  the  Kutta-Joukowski  law,  the  radial  com¬ 
ponents  of  vorticity  cause  both  a  normal  and  a  chordwise  loading  (in 
pounds  per  foot  of  rotor  radius).  However,  the  tangential  elements 
of  vorticity  and  the  radial  component  of  flow  velocity  cause  a  concen¬ 
trated  load  at  a  given  radial  position.  This  load  must  subsequently 
be  converted  into  a  distributed  loading  having  units  of  pounds  per  foot 
of  rotor  radius.  It  is  assumed  that  the  tangential  location  of  the  center 
of  pressure  for  each  radial  segment  is  the  same  as  that  measured  for 
the  two-dimensional  airfoil  section  at  the  angle  of  attack  of  that  segment. 
In  this  regard,  the  angle  of  attack  is  defined  in  the  plane  normal  to 
the  spanwise  axis  of  the  blade  segment.  Not  only  does  this  provide 
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a  valid  approximation  for  the  high-aspect-ratio  wing,  but  it  also 
yields  the  proper  tangential  location  for  the  center  of  pressure  in 
the  low-aspect-ratio  case. 

STRUCTURAL  MODEL 


Lagrangian  procedures  were  utilized  to  determine  the  equa¬ 
tions  of  motion  that  yield  the  accelerations  of  the  blades  in  their  rigid 
and  elastic  degrees  of  freedom.  The  resulting  inertial  and  elastic 
parameters  increase  geometrically  as  the  number  of  degrees  of 
freedom  used  to  represent  the  system  is  increased.  Accordingly, 
engineering  judgment  was  applied  to  reduce  the  number  of  terms 
that  result  from  a  mathematical  expansion  of  the  total  kinetic  and 
potential  energy  of  the  system.  This  reduction  was  accomplished  by 
dropping  terms  acknowledged  to  be  second  order  at  various  stages 
of  the  expansion.  However,  considerable  care  had  to  be  exercised 
in  order  not  to  oversimplify  prior  to  carrying  out  the  required  partial 
derivatives. 

An  excellent  example  of  the  magnitude  of  the  task  and  the 
complexity  of  the  result  that  is  obtained  upon  choosing  to  define  the 
kinetic  and  potential  energy  in  complete  detail  is  given  by  Lemnios 
in  Reference  [6],  Unfortunately,  the  equations  derived  in  this  latter 
work  could  not  be  reduced  to  the  equations  required  for  the  TBL-II 
study.  (In  particular,  the  equations  presented  in  Reference  [6] 
did  not  include  a  degree  of  freedom  corresponding  to  an  angular- 
displacement  response  of  the  rotor  shaft.  )  For  the  purposes  of 
this  study,  it  was  necessary  to  express  the  Lagrangian  of  the  system 


These  equations  are  applicable  to  rotors  assumed  to  have  a  flapping, 
feathering,  lagging,  and  servo-flap  deflection  degree  of  freedom, 
including  a  single  flatwise  bending  mode. 
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in  terms  of  four  rigid-blade  degrees  of  freedom  (i.  e.  ,  p  ,  fi  ,  <*,  and  9) 
and  a  finite  number  of  degrees  of  freedom  associated  with  blade 
flexibility.  On  assuming  that  a  flexible  blade  can  bend  in  both  the 
flatwise  and  chordwise  directions  and  can  also  twist  about  an  axis  defined 
as  the  elastic  axis,  a  choice  exists  for  representing  the  blades  by  means 
of  a  finite  number  of  lumped  mass  and  elastic  segments  or  by  means  of 
a  finite  number  of  assumed  radial  Junctions.  By  assuming  radial  functions 
that  correspond  to  the  resonant  vibration  modes  possessed  by  the  rotor 
blade  rotating  in  a  vacuum,  the  number  of  radial  functions  required  to 
represent  the  infinite  number  of  elastic  degrees  of  freedom  with  a  given 
level  of  accuracy  is  significantly  less  than  would  be  required  if  radial 
functions  were  assumed  on  a  completely  arbitrary  basis. 

For  purposes  of  the  TBL-1I  study,  the  elastic  response  of  a 
rotary  wing  was  assumed  to  be  described  by  a  finite  number  of  coupled 
flatwise-  and  chordwise-bending  modes  plus  a  finite  number  of  torsional 
modes.  (Although  the  number  of  assumed  modes  should  theoretically 
be  a  function  of  the  highest  frequency  present  in  the  disturbance  function, 
in  practice  this  number  is  very  often  cetermined  by  computational 
restrictions.  )  With  reference  tc  this  study,  it  is  clear  that  radial  func¬ 
tions  based  on  the  vibration  modes  exhibited  by  the  rotor  at  its  normal 
revolution  rate  should  give  accej/abie  answers  whenever  the  rotor  is 
rotating  at  this  rate  or  very  close  to  this  rate.  When  the  rotor  runs 
at  revolution  rates  other  than  the  normal  value,  the  resonant  frequencies 
and  associated  mode  shapes  arc  altered.  In  the  TBL-II  analysis,  it  is 
assumed  that  these  mode  shapes  Jo  not  change  to  an  appreciable  extent 
as  the  revolution  rate  of  the  rotor  is  varied.  Accordingly,  it  is  assumed 
that  the  mode  shapes  calculated  for  the  normal  revolution  rate  remain 
good  approximations  of  the  resonant  response  characteristics  of  the 
rotor.  Admittedly,  this  assumption  can  be  viewed  as  being  rather 
drastic  if  it  is  intended  to  make  transient-response  calculations  in 
which  the  rotor  approaches  or  achieves  a  near-zero  or  zero  revolution 
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rate.  It  shall  suffice  to  state  here  that  assumed  modes  of  constant 
character  are  used  to  represent  the  structural  dynamics  of  the 
rotor,  and  that  the  actual  extent  to  which  these  assumed  modes 
lead  to  inaccuracies  at  rotational  speeds  other  than  the  nominal 
value  cannot  be  stated  without  conducting  an  extensive  investigation 
of  this  aspect  of  the  problem.  Further,  it  can  be  argued  that  the 
assumption  of  mode  shapes  should  not  be  a  matter  of  concern  in  light 
of  the  assumptions  and  approximations  that  are  necessarily  made  to 
calculate  the  aerodynamic  forces  on  the  blade. 

In  deriving  the  Lsgrangian  of  the  rotor  system,  the  potential 
energy  of  the  system  is  assumed  to  derive  only  from  the  strain  energy 
existing  in  each  flexible  blade,  with  the  gravitational  energy  resulting 
from  the  displacements  of  the  rigid  blade  being  ignored.  Further,  no 
structural  damping  has  been  included  in  the  formulation  other  than  the 
structural  damping  assumed  to  exist  in  the  torsion  degree  of  freedom. 
It  should  be  noted  that  the  potential  energy  of  the  twisted-rotating 
beam  has  been  defined  in  accordance  with  the  theory  outlined  in 
Reference  [19],  as  further  discussed  in  Reference  [8].  In  view  of 
the  rather  complete  discussion  in  this  latter  reference,  the  details  of 
this  portion  of  the  analysis  are  not  included  in  this  report.  It  should 
be  noted,  however,  that  the  following  assumptions  have  been  made  in 
deriving  the  strain  energy  of  a  twisted,  rotating  beam: 

1.  Elementary  beam  theory  is  applicable. 

2.  Rotor-blade  cross  sections  are  symmetrical  with 
respect  to  the  major  neutral  axis;  hence,  the  c.  g. 
of  the  cross  section  lies  on  the  major  neutral  axis. 

3.  Gravity  forces  are  negligible  relative  to  centrifugal 
and  elastic  forces. 


Although  the  mode  shapes  remain  fixed,  the  modal  frequencies  are 
a  function  of  the  revolution  rate. 
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4.  Elastic  deformations  of  the  rotor  blade  are  small 


with  respect  to  the  radius  of  the  blade. 


The  derivation  of  the  equations  of  motion  for  an  articulated 
rotor,  which  has  the  hinge -offset  geometry  diagrammed  in  Figure  6, 
is  presented  in  Appendix  HI.  A  specification  of  the  boundary  conditions 
at  the  root  of  the  blade  and  a  designation  of  the  allowable  degrees  of 
freedom  cause  these  equations  to  reduce  to  a  form  applicable  to  the 
"rigid"  rotor  (i.  e. ,  no  fi  and  t?  degrees  of  freedom,  with  a  cantilever 
constraint  at  the  hub)  and  to  the  teetering  rotor  (i.  e. ,  no  degree  of 
freedom,  and  an  effective  cantilever  constraint  for  the  symmetric 
bending  modes  and  an  effective  pin  restraint  for  the  antisymmetric 
bending  modes,  including  the  flapping  mode). 

For  the  general  case  (i.e. ,  the  number  of  flexible  response 
modes  are  not  specified),  the  equations  of  Lagrange  yield  that: 
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generalized  force  corresponding  to  r** 
bending  mode  of  blade  i . 


=  generalized  force  corresponding  to  i*h 
torsion  mode  of  blade  i. 
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The  foregoing  equations  of  motion  for  blade  i  are  a  coupled 
set  since  the  inertial,  centrifugal,  and  elastic  forces  (as  expressed 
by  the  left-hand  side)  and  the  external  forces  (as  expressed  by  the  right- 

hand  side)  are  functions  of  all  of  the  degrees  of  freedom  possessed  by 
blade  i .  Further,  the  equations  for  blade  t  =  1  are  coupled  with  the 
equations  for  blades  i  =  2,  3,  4,  . .  .  through  the  if  degree  of  freedom 
which  is  common  to  all  blades.  Thus,  the  total  number  of  equations 
defining  the  response  of  an  articulated  rotor  is  equal  to  the  number  of 
rotor  blades  multiplied  by  the  assumed  number  of  blade  degrees  of 
freedom  plus  the  shaft  degree  of  freedom.  If,  for  example,  six  bending 
modes  and  two  torsion  modes  are  assumed  to  describe  the  elastic  re¬ 
sponse  of  a  blade,  then  an  articulated  blade  (i.  e.  ,  one  with  a  flapping, 
lagging,  and  feathering  degree  of  freedom)  would  possess  11  degrees 
of  freedom  plus  the  i//  degree  of  freedom.  A  four-bladed  rotor  would 
be  described  by  a  coupled  set  of  45  equations  of  motion;  a  two-bladed 
rotor  by  a  coupled  set  of  23  equations. 

In  order  to  minimize  the  numerical  difficulties  that  accompany 
the  simultaneous  solution  of  large  numbers  of  equations,  the  equation 
for  the  shaft-rotation  degree  of  freedom  was  reduced  to  its  simplest 
form  and  solved  independently.  The  resulting,  approximate  solution 
for  if  can  then  be  substituted  into  the  remaining  equations,  which  can 
be  solved  in  groups  corresponding  to  a  set  for  each  blade.  Thus, 
instead  of  determining  the  modal  accelerations  from  a  simultaneous 
solution  of,  say,  23  equations  for  a  two-bladed  rotor,  these  accel¬ 
erations  can  be  computed  by  performing  an  independent  solution 
oi  two  groups  of  equations  consisting  of  1 1  equations  per  blade. 

The  lead-lag  acceleration  calculated  for  each  blade  can  then  be  used 
to  correct  the  approximate  solution  for  the  angular  acceleration  of 
the  rotor  shaft, and  this  process  can  be  repeated  to  convergence. 

In  writing  the  equations  of  motion  for  a  rotor,  it  is  convenient 
to  use  a  matrix  notation  where  the  indices  m  and  i  indicate  an  instant 
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of  time  and  a  blade,  respectively, 
nate  a  generalized  coordinate,  i.  e. 


With  the  symbol  ^  used  to  desig- 


A  ‘ 

r 

C 

e 

. 

■  S  ■ 

u/  (t> 

. 

. 

<p"> 

i  j 

U(t)\ 

the  equations  of  motion  for  a  given  blade  can  be  represented  (symbolically) 

as 
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where 


[  •  L,  ■ 
{ L 


is  an  11  by  11  square  matrix  of  inertia  coefficients, 


is  a  column  matrix  of  11  generalized  accelerations. 


is  a  column  matrix  of  1 1  coefficients  defining  the 
inertial  coupling  between  |  jj  J  and  iy ,  and 


A 


is  a  column  matrix  of  11  quantities  representing 
the  combined  forces  resulting  from  centripetal 
and  Coriolis  accelerations,  elastic  deformations, 
and  aerodynamic  loadings. 


The  solution  of  Equation  (9),  yielding  generalized  accelerations,  veloc¬ 
ities, and  displacements,  can  be  expressed  as  follows: 
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Figure  12  is  a  block  diagram  of  the  manner  in  which  the  right- 
hand  side  of  Equation  (9)  was  subdivided  for  ease  of  computation.  The 
integrals  comprising  the  coefficients  in  the  indicated  square  and  column 
matrices  are  tabulated  in  Table  IV,  to  be  found  in  Appendix  III.  The 
external  force,  j  ,  shown  in  Figure  12,  is  the  generalized  force 

resulting  from  aerodynamic  loadings,  control -linkage  forces,  etc. 

The  integrals  yielding  the  generalized  forces  acting  on  each  blade  are 
given  in  Appendix  IV. 

On  assuming  six  modes  of  coupled  flatwise  and  chordwise 
bending,  and  two  modes  of  torsion  to  represent  the  elastic  response 
of  a  flapping,  lagging,  aid  feathering  blade,  one  finds  that  455  integrals 
must  be  evaluated  to  define  the  inertial,  centrifugal,  and  elastic  prop¬ 
erties  of  a  flexible,  articulated  blade.  The  tabulation  of  these  455 
integrals  has  been  omitted  from  this  report  since  they  represent  only 
a  detailed  expansion  of  the  integrals  defined  in  Table  IV. 

However,  it  is  pertinent  to  indicate  the  form  of  the  inertia 
matrix,  [>*]  ,  for  the  articulated  blade  prior  to  defining  the 
equations  of  motion  that  apply  to  a  two-bladed,  teetering  rotor.  If 
the  bending  and  flapping  modes  are  derived  as  orthogonal  modes,  and 
if  the  two  torsional  modes  are  assumed  to  be  orthogonal  to  each  other. 
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Figure  12.  TBL-H  EQUATIONS  OF  MOTION  IN  BLOCK-DIAGRAM  FORM. 


the  inertia  matrix  for  the  i  **  articulated  blade  has  the  form  indicated 
in  Figure  13.  Note  that  the  generalized  accelerations  ^  through 
are  bending  accelerations  and  the  accelerations  and 

are  torsional  accelerations. 


In  a  teetering  rotor,  the  lead-lag  degree  of  freedom  doe'?  not 
exist,  and  it  is  convenient  to  assume  that  the  j£ (¥),  <j(t* ,  and 
degrees  of  freedom  represent  antisymmetric  bending  modes  of  the 
pinned,  two-bladed  rotor.  On  letting  the  ,  and  degrees 

of  freedom  represent  symmetric  bending  modes  of  the  rotor,  and  on 
summing  the  inertial,  centrifugal,  and  elast'c  forces/moments  over 


the  entire  rotor  f~om  blade  tip  to  blade  tip,  equations  of  motion  for 
a  teetering  rotor  can  be  derived  in  terms  of  the  same  integrals  used 
for  the  articulated  blade.  These  equations  are  presented  in  Figure  14. 
It  should  be  noted  that  only  one  torsion  mode  has  been  assumed  to 
represent  the  elastic  twist  of  the  blades.  It  has  further  been  assumed 
that  each  blade  of  a  teetering  rotor  both  pitches  (against  a  flexible 
control  link)  and  twists  independently  of  the  other  blade.  With  these 
assumptions,  a  teetering  rotor  has  11  degrees  of  freedom,  as  is 
indicated  in  Figure  14.  The  kinematic  constraints  resulting  from 
the  assumed  division  of  symmetric  and  antisymmetric  modes  require 
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Figure  13  DEFINITION  OF  THE  INERTIA  MATRIX  FOR  THE  ith  ARTICULATED  BLADE. 
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Figure  14.  EQUATIONS  OF  MOTION  FOR  A  TEETERING  ROTOR. 


In  addition,  it  should  be  noted  that  t?  is  a  disturbance 
variable,and  therefore  the  initial  conditions  for  the  flap  angles,  y1/1  and 
,  must  reflect  a  preset  coning  angle,  if  such  an  angle  exists  in  a 
teetering  rotor. 
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COMPUTATIONAL  PROCEDURE 


GENERAL  REMARKS 

The  TBL-II  simulation  outlined  in  the  previous  section  was 
programmed  (with  the  aid  of  FORTRAN  IV)  for  solution  on  an  IBM-7044 
digital -computing  facility.  Th'  significant  expansion  in  the  number  of 
problem  variables,  relative  to  the  number  treated  in  TBL-I,  required 
that  calculations  be  restricted  to  a  two-bladed  rotor,  namely,  a  tee¬ 
tering  rotor  or  a  "rigid"  rotor.  In  general,  the  computational  method¬ 
ology  followed  that  employed  in  TBL-I.  Whereas  the  maximum  num¬ 
ber  of  radial  segments  used  to  represent  a  blade  in  TBL-I  was  10, 
TBL-II  provides  for  a  total  of  14  radial  segments.  This  increase 
was  deemed  necessary  to  obtain  accurate  loadings  at  the  root  of  a 
blade  during  high-advance-ratio  operations  in  which  the  root  of  the 
blade  becomes  the  leading  edge  of  a  surface  of  small  aspect  ratio. 

The  variable  angular  velocity,  if,  that  results  when  mechanical 
shaft  torque  is  an  input  quantity  integrates  to  nonrepeating  azimuth 
positions;  consequently,  it  is  more  difficult  to  judge  when  a  steady- 
state  operating  condition  is  achieved  than  was  the  case  for  TBL-I. 
Rather  than  employ  an  interpolation  scheme  to  determine  whether 
blade  response  variables  are  repeating  each  360  degrees  of  shaft 
rotation,  steady  state  was  defined  to  exist  when  the  mean  value  of 
(over  one  revolution  of  shaft  rotation,  approximate)  reaches  a  steady 
value. 


A  Runge-Kutta  integration  subroutine  is  employed  to  integrate 
the  computed  modal  accelerations  to  obtain  modal  velocities  and  dis¬ 
placements.  Selection  of  the  time  intervals  used  in  the  integration 
process  is  guided  by  a  rule  of  thumb  which  states  that  approximately 
10  time  intervals  are  required  per  cycle  of  response  in  order  to  obtain 
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an  accurate  integration.  The  time  interval  is  thus  governed  by  the 
natural  frequency  of  the  highest  frequency  mode  included  in  the 
simulation. 

In  order  to  compute  the  response  to  a  transient  input  of  shaft 
torque,  collective  pitch,  cyclic  pitch,  etc.  ,  it  is  necessary  to  integrate 
from  a  set  of  assumed  initial  conditions.  This  integration  must  proceed 
until  the  transient  response  resulting  from  the  disturbance  created  by 
the  assumed  initial  conditions  disappears  and  steady,  periodic  loadings 
and  blade  motions  are  achieved.  A  premium  is  thereby  placed  on  oper¬ 
ating  procedures  wherein  the  time  and  expense  of  achieving  an  initial 
steady  state  is  held  to  a  minimum. 

As  was  done  in  TBL-I,  the  program  is  first  operated  in  a 
mode  in  which  a  constant  induced  inflow  is  assumed.  With  the  wake 
model  inactive  (in  the  sense  that  induced-velocity  coefficients  are  not 
evaluated),  the  computing  time  required  to  integrate  from  an  assumed 
set  of  initial  conditions  to  an  approximate  steady-state  response  is 
significantly  reduced.  Upon  attainment  of  a  steady-state  response 
based  on  uniform  inflow,  the  wake  model  is  activated  and  computes 
induced  velocities  as  a  function  of  the  approximate  vortex  strengths 
that  were  previously  deposited  in  the  airstream.  The  integration 
process  continues  until  the  circulation  strengths  in  the  wake  become 
periodic  and  the  rotor  response  is  likewise  periodic. 

The  computing  time  required  for  one  calculation  interval  with 
the  wake  model  operative  c.nd  inoperative  are  9.  3  and  4.4  seconds 
respectively.  Although  this  saving  is  significant,  it  should  be  noted 
that  in  TBL-I  a  single  calculation  interval  required  only  2.  5  seconds 
of  machine  time  when  the  wake  model  was  inoperative.  This  difference 
between  TBL-I  and  TBL-II  is  attributable  to  the  larger  number  of  degrees 
of  freedom  and  integrations  associated  with  TBL-II. 
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Since  the  computing  time  associated  with  each  operation  of 
the  wake  model  is  quite  large,  the  computing  program  is  instructed  to 
calculate  a  wake  geometry  and  an  induced  velocity  in  the  plane  of  the 
rotor  at  specified  azimuth  intervals  rather  than  at  each  of  the  small 
time  intervals  required  to  perform  an  accurate  integration.  Provision 
is  also  made  for  specifying  a  delay  period  between  the  achievement  of 
a  steady  state  and  the  initiation  of  a  change  in  any  one  of  the  independent 
variables.  The  transient  variation  of  these  variables  (e.  g.  ,  shaft  torque 
or  collective  pitch)  is  treated  as  a  ramp  function  for  which  the  magnitude 
and  slope  can  be  specified. 

Various  aspects  of  the  computational  procedure,  other  than 
those  discussed  above,  are  the  same  as  were  employed  in  TBL-I  (see 
Reference  12). 

PROBLEMS  ENCOUNTERED 

Considerable  difficulties  were  encountered  in  starting  up  the 
developed  digital  simulation  when  all  of  the  degrees  of  freedom  discussed 
earlier  were  included  in  the  simulation.  Unfortunately,  insufficient  time 
was  available  to  determine  whether  the  resulting  unstable  responses  were 
a  result  of  programming  errors  or  were  caused  by  various  other  factors 
such  as  the  magnitude  of  the  disturbance  imposed  by  the  assumed  initial 
conditions.  Many  calculation  trials  were  made,  however,  and  these  have 
indicated  that  the  large  inertial  coupling  between  the  feathering  and 
torsion  degrees  of  freedom  places  a  very  severe  requirement  on  the 
initial  conditions  assumed  for  these  two  degrees  of  freedom.  If,  for 
example,  the  feathering  degree  of  freedom  is  assumed  to  possess  an 
initial  displacement  and  velocity  equivalent  to  that  of  the  swashplate, 
the  acceleration  computed  for  feathering  at  the  first  time  interval  is 
very  large,  resulting  in  a  torsional  moment  (from  the  large  inertial 
coupling)  that,  in  turn,  causes  a  large  torsional  acceleration.  If  the 
initial  departure  from  dynamic  equilibrium  is  large,  as  was  evidenced 
by  the  feathering  and  torsional  accelerations  computed  at  the  first  time 
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interval,  the  rotor  response  proved  to  be  unstable  as  a  result  of  the 
direct  influence  of  feathering  and  torsion  on  blade  angles  of  attack  and 
on  the  associated  aerodynamic  loading. 

The  expedient  adopted  to  alleviate  this  difficulty  consisted  of  the 
effective  removal  of  the  compliant  control  link  between  the  swashplate 
and  the  blade,  by  this  means,  the  feathering  degree  of  freedom  becomes 
an  independent  input  variable  rather  than  a  response  variable.  Since  the 
feathering  accelerations  are  then  prescribed  to  be  the  values  determined 
by  the  tilt  of  the  swashplate,  the  inertial  coupling  between  torsion  and 
feathering  is  automatically  bounded  to  reasonable  values.  The  dynamic 
systf  .1  obtained  on  removal  of  the  feathering  degree  of  freedom  exhibited 
less  sensitivity  to  assumed  initial  conditions  and,  accordingly,  was  em- 
ployed  in  the  remainder  of  this  study. 

Another  computational  matter  worth  discussing  at  this  point 
relates  to  the  determination  of  the  unknown  circulations  from  the  simul¬ 
taneous  equations  defining  the  boundary  condition  of  flow  tangency.  Prior 
to  the  development  of  the  aerodynamic  model  based  on  the  abandonment 
of  the  Prandtl  simplification,  the  solution  of  the  simultaneous  equations 
was  iterated  until  the  boundary  condition  was  satisfied  at  all  control 
points,  even  though  some  control  points  were  located  on  blade  segments 
experiencing  angles  of  attack  in  excess  of  stall.  In  the  calculations  that 
utilized  two-dimensional  section  data,  this  procedure  was  implemented 
by  taking  the  solution  of  the  linear  equations,  computing  the  resulting 


A  subsequent  discovery  of  the  crucial  importance  of  the  quasi-static 
aerodynamic  moment  (caused  by  pitching  and  torsional  velocities)  in 
providing  adequate  damping  for  the  stability  of  the  torsional  mode  indi 
cates  that  the  above -discussed  problem  may  not  have  existed  had  this 
damping  moment  been  included  in  the  equations  of  motion  from  the 
beginning  of  the  study. 
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angles  of  attack,  redefining  the  circulations,  recomputing  the  angles 
of  attack,  etc.,  until  convergence  was  obtained. 

In  setting  up  the  current  procedure  for  computing  the  unknown 
circulations,  it  was  recognized  that  there  is  an  inconsistency  in  demanding 
that  there  be  flow  tangency  at  those  sections  where  the  blade  is  presum¬ 
ably  stalled.  Accordingly,  a  procedure  was  adopted  in  which  the  linear 
equations  were  solved  once  and  the  resulting  circulations  were  truncated 
if  the  local  two-dimensional  angle  of  attack  exceeded  a  stall  ang.'e  computed 
as  a  function  of  the  local  Mach  number.  The  local  two-dimensional  angle 
of  attack  was  first  determined  as  a  function  of  (a)  blade  geometry  and 
motion  and  (b)  the  velocity  induced  by  the  "old"  wake.  The  additional 
velocity  induced  by  the  adjacent  wake  elements  was  then  computed  using 
the  truncated  circulations  and,  after  making  a  single  correction  to  the 
local  Mach  number  and  the  local  angle  of  attack,  the  truncation  piocess 
was  repeated.  No  attempt  was  made  to  iterate  to  values  of  circulation 
that  would  insure  flow  tangency  everywhere.  Lacking  any  theory  or 
rationale  which  indicates  how  the  circulation  on  the  blades  and  in  the 
adjacent  wake  elements  should  readjust  their  strengths  whenever  certain 
spanwise  segments  are  stalled,  there  did  not  appear  to  be  any  alternative 
other  than  to  accept  the  above-described  "open-loop"  solution  procedure. 
This  decision  did,  however,  produce  an  unanticipated  result;  namely, 
a  fluctuation  in  calculated  loadings  where  the  fluctuation  is  caused  by 
the  method  of  solution  and  results  intermittently  in  oscillatory  loadings, 
with  the  frequency  of  this  oscillation  being  the  calculation  frequency 
(one  calculation  every  0.004  second  or  250  calculation  intervals  per 
second).  Fortunately,  this  frequency  is  sufficiently  high  that  the  conse¬ 
quences  of  this  spurious  high-frequency  forcing  are  negligible, in  that 
this  high-frequency  excitation  does  not  produce  a  significant  response 
in  any  of  the  assumed  modes. 
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INPUT  DATA 

A  separate  computational  program  was  used  to  evaluate  the 
numerous  integrals  defined  in  the  TBL-II  analysis.  Specific  input 
data  were  required,  making  it  desirable  to  select  a  rotor  for  which 
the  necessary  structural  information  was  available.  The  further 
restrictions  imposed  by  the  capacity  of  available  computing  equip¬ 
ment  led  to  the  selection  of  the  two-bladed  rotor  of  the  HU-1A 
helicopter. 
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Blade  cross-section  geometry  and  material  properties,  as 
supplied  by  the  Bell  Helicopter  Company  (Reference  [22]),  were 
analyzed  at  1 1  spanwise  stations  to  obtain  the  structural  parameters 
defined  in  Appendix  III.  These  results  are  tabulated  in  Table  I. 

Existing  data  on  the  mass,  *>($),  and  bending -stiffness,  £TXX  and 

distributions  (see  Figure  15)  were  employed  in  a  modal  analysis 
(Reference  [20])  to  obtain  mode  shapes  for  three-symmetric  and  three- 
antisymmetric  bending  modes.  The  coupled  flapwise-chordwise  bending 
mode  shapes  obtained  for  the  nominal  rotational  speed  of  314  rpm  are 
presented  in  Figure  16,  together  with  the  two  uncoupled  torsional-mode 
shapes  calculated  by  means  of  a  Holzer  analysis.  Figure  17  shows  the 
straight  elastic  axis  that  was  assumed  to  approximate  the  computed 
elastic -axis  location  (see  Table  I)  and  also  shows  the  location  of  the 
local  c.  g.  ,  the  one-quarter  chord,  and  the  feathering  axis  relative  to 
the  elastic  axis. 


A  lack  of  familiarity  with  the  vibratory  behavior  of  highly 
twisted,  rotating  blades,  as  yielded  by  a  coupled  flatwise-  and  chord- 
wise-bending  analysis,  provided  cause  for  questioning  the  data  plotted 
in  Figure  16.  In  order  to  produce  a  higher  level  of  confidence  in 
these  results,  calculations  were  made  for  a  range  of  rotational  speeds. 
These  calculations  served  to  demonstrate  the  manner  in  which  the 


48 


TABLE  I 

STRUCTURAL  PARAMETERS  OF  THE  HU-1  A 
BLADE  DETERMINED  FROM  CROSS-SECTION  DATA 


SPAHNISE 

STATION 

in. 

A 

ft2 

ft 

ft2 

ft6 

ft* 

M- 

a  lug* -ft2 
ft 

42 

5.62  x  10’2 

0.472 

0.133 

-2.07  x  I0‘3 

1.31  x  I0"3 

0.04008 

55 

4.14  x  10'2 

0.451 

0.155 

-1.99  x  10'® 

1.17  x  10”3 

0.03440 

65 

3.50  x  10"2 

0.441 

0.165 

-1.86  x  10"3 

1.03  x  10-3 

0.03109 

72.5 

3.26  x  I0‘2 

0.443 

0.165 

-1.66  x  10’3 

0.946x  10-3 

0.02874 

65 

3.31  x  I0-2 

0.440 

0.163 

-1.79  x  )0‘3 

1.01  x  10"1 

0.02687 

IOO 

3.10  x  10"2 

0.412 

0.146 

-1.67  x  10"® 

1.00  x  I0’3 

0.02541 

120 

3.06  x  I0”2 

0.381 

0.132 

-1.41  x  10-3 

0.839x  tO-3 

0.02160 

140 

2.95  x  10'2 

0.364 

0.116 

-1.20  x  10‘3 

0.796x  10”3 

0.01640 

•  '65 

2.97  x  I0"2 

0.347 

0.116 

-1.29  x  10’3 

0.649x  10-3 

0.01913 

170 

3.16  x  I0’2 

0.336 

0.114 

-1.39  x  10"3 

0.870x  I0“3 

0.02346 

240 

4.31  x  I0-2 

0.270 

0.102 

-2.06  x  10"3 

1.45  x  10"3 

0.02156 

NOTE:  SEE  APPENDIX  III  FOR  THE  DEFINITION  OF  4,  «4,  k\ ,  3„  AND 
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OCATION  OF  ELASTIC  AXIS,  SECTIONAL  CENTER  OF  GRAVITY,  AND  FEATHERING 
tXIS  WITH  RESPECT  TO  THE  LEADING  EDGE  OF  THE  HU-14  8LADE. 


coupled  bending  modes  chat  ge  their  character  as  rotational  speed 
is  varied  and  provided  assurance  that  the  mode  shapes  computed 
for  the  314-rpm  case  were  realistic. 

The  modal  frequencies  obtained  in  these  calculations  are 
plotted  in  Figure  18.  Note  that  both  the  first  and  second  symmetric 
and  antisymmetric  bending  modes  change  from  a  flatwise  and  chord - 
wise  mode  to  chordwise  and  flatwise  modes,  respectively.  It  is  seen 
that  the  first  and  second  symmetric  modes  approach  each  other  in 
frequency  and  also  change  their  chordwise  to  flatwise  displacement 
ratio  in  the  vicinity  of  314  rpm,  the  nominal  rotational  speed  of  the 
HU-iA  rotor.  This  particular  result  does  raise  a  question  as  to 
whether  uncoupled  bending  modes  might  have  been  better  choices 
for  assumed  modes,  since  it  appears,  in  theory,  that  it  is  preferable 
to  assume  mode  shapes  that  do  not  change  appreciably  as  the  rotational 
speed  of  the  rotor  is  varied.  The  data  provided  by  the  present  program 
are  not  sufficient  to  determine  a  preferred  procedure  for  choosing 
"assumed  modes". 


ROTOR,  rpa 

RESONANT  FREQUENCIES  OF  COUPLED  FLATWISE  AND  CHORDWISE 
BENDING  VERSUS  ROTOR  RPM;  HU-1  A  ROTOR  BLADES. 

(ROOT  CHORD  IS  AT  ZER0-DE6REE  COLLECTIVE  PITCH) 


EQUIVALENCE  BETWEEN  THE  TWO  AERODYNAMIC  MODELS  AT 
LOW  ADVANCE  RATIOS 


In  the  section  dealing  with  the  aerodynamic  model,  it  was 
shown  that  the  vortex  geometry,  adopted  in  TBL-H  to  represent  the 
bound  vorticity  on  the  blade,  yields  a  realistic  distribution  of  loadings 
on  a  low-aspect-ratio  surface  corresponding  to  a  two-bladed  rotor 
stopped  at  the  0  -  180-degree  azimuth  position.  As  in  discussed  in 
Appendix  II,  features  were  included  in  the  model  in  order  to  yield,  at 
low  advance  ratios,  circulation  strengths  equivalent  to  the  lift  coefficient 
that  is  obtained  from  a  "strip"  analysis  employing  two -dimensional 
empirical  airfoil  data.  In  view  of  the  heuristic  nature  of  the  steps  taken 
to  insure  matching  results  from  the  two  models  at  low  advance  ratios, 
questions  did  exist  as  to  the  level  of  agreement  that  would  be  obtained  in 
actual  computational  practice. 

Figure  19  shows  that,  for  a r>  advancing  blade  (/t  =0.26)  operating 
\t  angles  of  attack  below  the  stall,  excellent  agreement  is  obtained  be¬ 
tween  the  spanwise  distribution  of  loadings  yielded  by  the  two  aerodynamic 
models.  The  minor  differences  noted  at  the  root  of  the  blade  are  due,  in 
part,  to  the  assumption  that  the  negative  component  of  radial  flow  expe¬ 
rienced  at  f  =142  degrees  causes  the  root  to  be,  in  effect,  a  trailing 
edge.  The  retreating  blade,  i.  e. ,  =  322  degrees,  happens  to  be 

encountering  angles  of  attack  in  excess  of  the  stall  angle;  therefore, 
the  loadings  computed  on  the  basis  of  two-dimensional  section  data  are 
less  than  those  obtained  from  the  solutions  of  the  equations  expressing 
the  boundary  condition,  where  these  solutions  have  been  limited  to  a 
specified  value  of  maximum  circulation.  The  differences  noted  for  the 
retreating  blade  are,  therefore,  to  be  expected.  It  was  concluded  that 
the  agreement  between  the  two  aerodynamic  models  is  quite  adequate  at 
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Figure  19 


COMPARISON  OF  SPANWI SE  DISTRIBUTION  OF  CIRCULATION  YIELDED 
BY  TWO  AERODYNAMIC  MODELS,  u.  =  0.26. 
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low  advance  ratios  and,  therefore,  it  is  postulated  that  the  aerodynamic 
analysis  employed  in  this  study  has  a  rational  basis  for  predicting  load 
distributions  on  long,  narrow  blades  irrespective  of  the  local  orientation 
of  the  total  velocity  vector  and  the  magnitude  of  the  radial  component  of 
this  velocity  vector. 

Calculations  made  for  an  assumed  rigid  rotor  stopped  at  the 
0  -  180- degree  azimuth  position  have  pointed  up,  however,  certain 
shortcomings  of  the  wake  model  employed  in  the  study.  In  particular, 
the  methods  used  to  define  the  location  and  strength  of  the  rolled-up  tip 
and  root  vortices  produced  at  low  advance  ratios  do  not  carry  over,  in 
a  logical  manner,  for  representing  the  roll-up  process  behind  a  low- 
aspect-ratio  surface.  The  wake  trailed  behind  the  high -aspect-ratio 
surface  is  a  function  of  the  change  in  circulation  (integrated  over  the 
chord)  along  the  s^an,  with  the  strengths  of  the  root  and  tip  vortex  being 
equivalent  to  the  maximum  circulation  achieved  along  the  span.  For  the 
low-aspect-ratio  surface,  the  adopted  horseshoe-vortex  representation 
yields  the  distribution  of  vorticity  along  the  "chord"  dimension.  Con¬ 
sequently,  the  logic  for  assigning  strengths  to  rolled-up  elements  of  the 
sheet  of  vorticity  trailed  from  the  side  edges  of  the  airfoil  is  incorrect. 

Note  that  the  variation  in  the  strength  of  this  vortex  sheet  is  equivalent 
to  the  chordwise  distribution  of  bound  vorticity  and  that  the  strength  of 
a  single  vortex  element,  representing  a  roll-up  of  this  sheet,  should  be 
the  integral  of  the  chordwise  distribution  of  vorticity.  This  is  not  the 
case  in  the  adopted  model.  It  appears  that  although  the  adjacent  non- 
rolle<'  elements  of  wake  streaming  back  from  the  side  edge  of  a  low-aspect- 
ratio  surface  yields  an  appropriate  geometry  for  representing  its  chord- 
wise  vorticity  distribution,  the  rolled  elements  (two  in  number)  have  little 
justification,  either  with  respect  to  their  location  or  their  strengths. 
Fortunately,  the  logic  which  is  inappropriate  to  the  low-aspect-ratio 
flow  case  gives  these  rolled-up  elements  a  geometry  and  strength  which, 
for  the  most  part,  make  their  influence  self-canceling.  In  order  to 
eliminate  the  steady-state  loading  error  caused  by  this  self-imposed 
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artifact,  the  steady-state  computation  procedure  was  modified  to  effec¬ 
tively  eliminate  the  rolling-up  process  while  the  rotor  was  stopped  at 
the  0  -  183-degree  azimuth  position.  On  initiating  the  transient,  that 
is,  engine  torque  is  applied  to  start  up  the  rotor,  the  rolling -up 
process  necessarily  reenters  the  calculation  procedure.  The  wake 
geometry  resulting  from  the  roll-up  logic  built  into  the  calculation  pro¬ 
cedure  accordingly  contributes  an  unknown  level  of  error  in  calculating 
loadings  at  very  high  advance  ratios  when  the  rotor  is  in  the  vicinity  of 
0-180  degrees  of  azimuth.  At  azimuth  locations  in  the  vicinity  of  90  - 
270  degrees,  the  wake  model  operates  in  a  satisfactory  manner,  irre¬ 
spective  of  the  magnitude  of  the  advance  ratio. 

RESPONSE  TO  POWER  LOSS  AND  A  BRAKING  TORQUE 

Three  transient  responses  were  computed  to  examine  the  aero¬ 
dynamic  loadings  and  the  structural  motions  of  a  two-bladed  teetering 
and  rigid  rotor  experiencing  a  power  loss  at  a  high-speed  cruise  condition. 
In  order  to  obtain  representative  input  data  for  shaft  angle,  cyclic  pitch, 
etc.,  experimental  values  reported  in  Reference  [24]  for  a  forward  flight 
speed  of  111.3  knots  (/<.  =  0.261)  were  employed.  The  tip-path  plane 
computed  for  the  teetering-rotor  configuration  (based  on  blade  flapping 
only)  was  used  to  determine  a  shaft  angle  which  would  give  the  rigid  rotor 
the  same  tip -path-plane  orientation.  In  adjusting  collective  pitch  on  the 
rigid  rotor  to  yield  the  same  average  lift  produced  by  the  teetering  rotor, 
slight  differences  in  average  lift  were  accepted  resulting  in  different 
steady-state  torque  requirements  as  tabulated  in  Table  II. 

Transient  responses  caused  by  both  a  50-  and  100-percent  power 
loss  were  computed  for  the  teetering-rotor  configuration.  Figure  20  indi¬ 
cates  that  the  50-percent  power  loss  was  applied  in  0.  1  second;  that  is, 
in  approximately  one -half  of  a  rotor  revolution.  The  time  history  in 
Figure  20  is  limited  to  1.  3  seconds  and  shows  one  complete  revolution 
of  steady-state  behavior  prior  to  initiating  the  reduction  in  shaft  torque. 
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RESPONSE  OF  TEETERING  ROTOR  TO  SO%  POWER  LOSS;  ^  =  0.26,  9e  =  0.369  RADIAN,  as  =  0.113  RADIAN 
(a)  AERODYNAMIC  LOAD  AND  FORCE  TIME  Hi  STORIES. 


RESPONSE  OF  TEETERING  ROTOR  TO  50*  POWER  LOSS;  /*  =  0.26,  9e  =  0-369  RADIAN,  <3*.=  0.113  RADIAN 
(b)  BLADE  DISPLACEMENTS  AND  BENDING  MOMENTS;  SHAFT-VELOCITY  RESPONSE. 


TABLE  II 

STEADY-STATE  TRIM  CONDITIONS; 
POWER  LOSS  AND  BRAKING  TRANSIENTS 


ROTOR 

CONFIGURATION 

TEETERING 

(FULL  THRUST) 

RIGID 

(FULL  THRUST) 

RIGID 

(ZERO  THRUST' 

V 

-  fp* 

188 

188 

168 

as 

-  rsd 

0. 1 1395 

0.10400 

0 

&c 

-  rsd 

0.360 

0.352 

0.138 

o„ 

-  rad 

-0.1278 

-0.1278 

-0.00672 

0,C 

-  rad 

0.0408 

0.0408 

0.00417 

a 

-  rp« 

313 

313 

313 

M 

0.261 

0.261 

0.260 

Qo 

-  ft-lb 

7,600 

6,650 

1,800 

The  oscillations  in  loading  occurring  at  the  calculation  frequency  are 
spurious,  as  has  been  discussed  earlier.  In  most  cases,  it  is  possible 

to  distinguish  between  loading  "spikes"  that  derive  from  sudden  changes 
in  the  induced-velocity  field  and  the  oscillations  that  derive  from  the 
calculation  procedure.  In  general,  all  of  the  plotted  responses  are  self- 
explanatory.  However,  it  should  be  noted  that  "Aerodynamic  Shear- Y 
and  Z  Component"  refers  to  the  hub  shear  produced  by  aerodynamic 
loadings  on  one  blade  as  measured  in  the  shaft  plane  and  perpendicular 
to  the  shaft  plane,  respectively. 

Examination  of  the  loadings  plotted  for  five  span  locations  does 
not  reveal  any  unusual  or  surprising  behavior.  It  is  seen,  however,  that 
the  mean  response  for  the  symmetric  bending  modes  is  much  larger  than 
would  normally  be  anticipated  for  a  preconed  rotor.  Investigation  has 
shown  that  the  conversion  of  the  equations  applicable  to  an  articulated  rotor 
to  those  applicable  to  a  teetering  rotor  has  inadvertently  caused  the  influ¬ 
ence  of  preconing  on  symmetric  bending  to  be  omitted.  Specifically,  it 
was  found  that  the  assumed  orthogonality  between  flapping  and  symmetric 
bending  does  not  hold  for  the  steady,  preconed  component  of  flap  angle. 
Although  this  omission  of  the  influence  of  preconing  on  steady  bending  was 
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not  discovered  until  all  the  calculations  had  been  completed,  it  is  be¬ 
lieved  that  conclusions  drawn  with  respect  to  either  dynamic  and/or 
transient  phenomena  are  not  significantly  altered  by  the  large  values 
of  mean  bending  and  mean-bending  moments  that  have  been  computed. 

This  comment  is  believed  to  hold  for  the  other  cases  studied  as  well 
as  for  the  response  time  histories  under  discussion.  In  effect,  the 
structural  responses  computed  in  this  investigation  apply  to  an  unconed 
teetering  rotor;  the  aerodynamic  forces  do,  however,  reflect  the  in¬ 
fluence  of  3  degrees  of  coning  angle. 

Figure  20  shows  that  the  peak-to-peak  amplitude  of  the  sym¬ 
metric-bending  modes  do  not  change  appreciably  following  the  power 
loss, although  the  mean  deflections  of  the  symmetric  modes  decrease 
with  the  decreased  loading  as  would  be  expected.  The  peak-to-peak 
response  in  flappirg  and  in  the  first  antisymmetric  mode  increases 
with  time  as  the  rotor  begins  to  slow  down.  The  decrease  in  overall 
loadings  causes  a  decrease  in  the  mean  flatwise  bending  moment  (at 
10-percent  span)  whereas  the  torque  disturbance  causes  a  significant 
increase  in  the  peak-to-peak  chordwise  deflection  of  the  blades  re¬ 
sulting  in  an  equivalent  increase  in  the  peak-to-peak  chordwise  bending 
moments.  It  should  be  noted  that  loading  disturbances  which  cause 
increased  excitation  of  the  torsion  mode  (25.  8  cps)  result  also  in 
increased  excitation  of  the  third  antisymmetric  bending  mode  (25.  1  cps). 

Figure  21  shows  that  the  100 -percent  power  loss  was  initiated 
at  an  azimuth  position  of  270  degrees, whereas  the  50 -percent  power  loss 
occurred  at  120  degrees  of  azimuth.  Further,  the  power  loss  was  applied 
at  the  same  rate  as  previously,  thus  requiring  0.  2  second,  or  approxi¬ 
mately  one  rotor  revolution,  to  lose  full  power.  The  response  produced 
by  the  100-percent  power  loss  is  similar  to  that  produced  by  the  50- 
percent  power  loss.  However,  it  appears  that  greater  high-frequency 
torsional  and  high-frequency  flatwise-bending-moment  responses  are 
excited  by  the  100 -percent  power  loss.  In  addition,  the  time  history  is 
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RESPONSE  OF  TEETERING  ROTOR  TO  100%  POWER  LOSS;,*  =  0-26,  =  0.369  RADIAN,  as  =  0,113  RADIAN 

(a)  AERODYNAMIC  LOAD  AND  FORCE  TIME  HISTORIES. 
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RESPONSE  OF  TEETERING  ROTOR  TO  100^  POWER  LOSS;^  =  0.26,  Bc  =  0  .  369  R/ 
(b)  BLADE  DISPLACEMENTS  AND  BENDING  MOMENTS;  SHAFT-VELOC I TY  RESPONSE 


long  enough  to  show  that  the  oscillation  in  the  chordwise  bending  moment 
(excited  by  the  change  in  shaft  torque)  finally  begins  to  decay  about 
0.  9  second  after  completion  of  the  input  disturbance. 


Figure  22  indicates  that  loadings  on  the  rigid  rotor  are  very 
similar  to  those  produced  on  the  teetering  rotor.  The  bending  and 
torsional  responses  of  the  rigid  blades  are,  however,  significantly 
different  from  those  experienced  by  the  teetering  rotor.  Whereas 
first  and  second  symmetric  bendings  on  the  teetering  rotor  are  in 
phase  with  each  other,  resulting  in  increased  flatwise  displacements 
accompanied  by  decreased  chordwise  displacements,  the  first  and 
second  cantilever  modes  on  the  rigid  rotor  tend  to  be  out  of  phase, 
with  smaller  components  of  second  harmonic  response  being  present. 
Mothwithstanding  the  out-of-phase  addition  of  first  and  second  bending, 
the  resulting  peak-to-peak  excursions  in  flatwise  and  chordwise  bending 
displacements  and  bending  moments  are  substantially  greater  than  those 
experienced  by  the  teetering  rotor.  In  fact,  the  periodic  chordwise 
bending  moment  for  the  rigid  rotor  during  steady-state  flight  is  signif¬ 
icantly  greater  than  the  chordwise  bending  -  moment  transient  resulting 
from  a  power  loss  on  the  teetering  rotor. 

It  should  be  noted  that  Figures  20,  21,  and  22  present  only  a 
small  portion  of  the  data  that  are  or  can  be  obtained  in  the  transient- 
response  calculations  performed  in  this  study.  Although  the  original 
plan  called  for  presenting  the  inertial  hub  shears  produced  by  each 
blade  in  addition  to  the  aerodynamic  shears,  it  was  discovered  that  a 
required  axis  transformation  was  not  introduced  into  the  computation 
of  the  inertial  shears;  consequently,  the  computed  results  were  in 
error.  (This  discovery  was  made  sufficiently  late  in  the  program  that 
time  did  not  permit  a  correction  of  the  noted  deficiency.  )  Although  the 
presented  time  histories  give  an  adequate  accounting  of  the  loading 
changes  that  occur  as  a  function  of  the  azimuthal  position  of  the  blade, 
the  spanwise  distribution  of  loadings  and  the  corresponding  spanwise 


69 


distribution  of  bending  displacements  and  moments  are  not  adequately 
documented  in  this  type  of  presentation.  The  spanwise  distribution 
of  loadings  and  the  bent  shape  of  the  blade  (flatwise  and  chordwise) 
are  available  in  printed  tabulations,  but  their  bulk  precludes  their 
presentation  in  this  report.  It  should  be  pointed  out  that  the  flatwise 
and  chordwise  bending  moments  were  computed  at  only  one  span 
position,  namely,  10  percent.  It  is  felt  that  a  complete  examination 
of  the  influence  of  power  loss  and  braking  torques,  in  terms  of  the 
transient  stresses  imposed  on  rotor  blades,  would  require,  moreover, 
a  thorough  review  of  the  spacial  distribution  of  bending  moments  (flat¬ 
wise  and  chordwise). 

The  response  to  a  rapid  and  severe  braking  torque  applied  to 
the  shaft  of  a  rotor  (operating  at  the  zero-lift  condition)  is  presented  in 
Figure  23.  In  this  instance,  the  wake  is  streaming  back  from  the  blades, 
remaining  essentially  in  the  plane  of  the  rotor.  The  operation  condition 
(defined  in  Table  II)  is  assumed  to  be  typical  of  a  rotor  mounted  on  a 
compound  helicopter  which  is  being  supported  by  the  lift  on  its  fixed  wing, 
just  prior  to  removing  power  (required  to  turn  the  rotor  at  100-percent 
rpm),  and  then  braking  the  rotor  to  a  stop.  The  impulses  in  loading 
caused  by  the  close  proximity  of  the  wake  to  the  blades  are  quite  evident. 
It  shou1J  be  noted  that  some  of  these  impulses  are  quite  severe  and 
constitute  a  very  irregular  spanwise  distribution  of  loading. 

When  the  collective  pitch  is  set  to  yield  a  mean,  nominal  value 
of  zero  thrust,  the  loading  outboard  on  the  blades  is  negative,  with  the 
loading  inboard  being  positive.  As  the  rotor  slows  down,  the  negative 
loads  outboard  are  reduced  and  the  total  rotor  lift  is  seen  to  increase  to 
a  small,  mean,  positive  value.  At  the  end  of  the  time  history  shown  in 
the  figure,  the  angular  velocity  of  the  rotor  has  dropped  from  an  initial 
value  of  32.81  radians  per  second  to  18.  98  radians  per  second.  Thus, 
the  advance  ratio,  ,  has  increased  from  0.  26  to  0.  45.  Within  this 
range  of  advance  ratio,  it  appears  that  no  unusual  loading  phenomena 
have  been  encountered  other  than  the  wake-disturbance  effects  that 
result  from  the  absence  of  the  usual  downward  transport  velocity. 
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Figure  22.  RESPONSE  OF  RIGID  ROTOR  TO  100%  POWER  LOSS;  ^  =  0.26,  Gc  =  0.352  RADIAN,  <*s  =  0.101  RADIAN 
(a)  AERODYNAMIC  LOAD  AND  FORCE  TIME  HISTORIES. 
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Figure  22.  RESPONSE  OF  RIGID  ROTOR  TO  100%  POWER  LOSS;  =  0.26,  6C  =  0.352  RADIAN,  <^=  0.104  RADIAN 
(b)  BLADE  DISPLACEMENTS  AND  BENDING  MOMENTS;  SHAFT- VELOCITY  RESPONSE. 
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Figure  23.  RESPONSE  OF  RIGID  ROTOR  AT  ZERO  THRUST  TO  LARGE  BRAKING  TORQUE; 

INITIAL  ADVANCE  RATIO  =  0.26,  9C  =  0.138  RADIAN,  as=  0  RADI Ari. 

(b)  BLADE  01 Sf  .ACEMENTS  AND  BENDING  MOMENTS;  SHAFT- VELOCITY  RESPONSE. 


LOADINGS  AND  BLADE  RESPONSE  EXPERIENCED  BY  A  RIGID  ROTOR 
IN  THE  "STOPPED"  AND  "START-UP”  CONDITIONS 

Figure  24  presents  the  loadings  and  displacements  experienced 
by  the  rotor  blade  when  it  is  stopped  in  the  downstream  position  (0-degr*.e 
azimuth)  and  then  slowly  rotated  through  one  full  rotor  revolution.  In  this 
instance,  it  would  have  been  desirable  to  plot  also  the  results  obtained  for 
the  blade  which  is  stopped  in  the  upstream  position.  However,  the  tape 
provided  to  the  automatic -plotting  routine  contained  results  only  for  the 
reference  blade;  that  is,  the  blade  which  is  located  at  f  =  0  degree  when 
the  rotor  is  assumed  to  be  at  <p  -  0  degree. 

It  should  be  realized  that  the  results  shown  for  i  =  0  second  have 
been  obtained  by  assuming  zero  values  of  blade  displacement  and  velocity 
for  the  initial  conditions.  Since  gravitational  forces  have  been  omitted, 
the  blade  does  not  bend  due  to  its  own  weight  and  this  means  that  the  blade 
is  initially  straight  when  first  exposed  to  the  airstream.  Since  the  as¬ 
sumed  rigid  rotor  has  3  degrees  of  preconing  and  the  shaft  angle  of 
attack  is  assumed  to  be  zero,  the  forward  blade  is  at  3  degrees  positive 
angle  of  attack  with  respect  to  the  airstream,  and  the  downstream 
blade  has  a  3 -degree  negative  angle  of  attack  with  respect  to  the 
airstream.  The  loads  and  blade  bending  resulting  from  these  initial 
conditions  are  permitted  to  reach  a  steady,  equilibrium  state  prior  to 
"unstopping"  the  rotor  and  applying  a  positive  shaft  torque  to  speed  up 
the  rotor. 


!  . 
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Examination  of  Figure  24  indicates  that  any  prior  fears  of  the 
possibility  of  static  divergence  or  dynamic  instability  appear  to  be 
unwarranted.  Obviously,  there  is  need  to  proceed  to  higher  rotational 


It  should  be  noted,  however,  that  these  blades  have  12  degrees  of 
negative  twist  with  the  collective  pitch  (as  measured  at  the  root  chord) 
set  to  approximately  7.  9  degrees. 
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speeds  and  to  consider  the  changes  that  may  result  from  representing 
the  mode  shapes  more  accurately  than  has  been  done  here,  before 
drawing  firm  and  final  conclusions  on  this  point.  Granting  that  the 
approximate  aerodynamic  analysis  employed  leads  to  inexact  definition 
of  the  loadings,  the  evidence  appears  to  support  the  conclusion  that  this 
specific  rotor  configuration  is  able  to  stop  and  start  up  without  encoun¬ 
tering  an  overstress  condition. 

It  is  instructive  to  examine  Figure  24  in  greater  detail.  Although, 
at  *fJ  =  0  degree,  the  wake  from  the  upstream  blade  is  passing  over  the 
blade  whose  loading  is  plotted  in  this  figure,  it  is  seen  that,  notwith¬ 
standing  this  disturbance  factor,  the  loading  has  the  expected  negative 
values  at  the  2.  5-  and  12.  5-percent  span  locations  which  correspond  to  the 
forward  chordwise  portions  of  the  resulting  low-aspect-ratio  surface.  It 
appears  that  in  the  neighborhood  of  ^  =  5  degrees,  the  control  point  of 
the  87.  5-percent  span  station  passes  in  close  proximity  to  a  wake  element 
originating  from  the  forward  blade  as  a  result  of  the  wake -definition 
logic  which  places  a  rolled-up  element  in  space  as  if  it  had  originated 
from  the  trailing  edge  (90-percent  span  location),  thereby  causing  this 
rolled-up  element  to  be  located  to  the  right  of  the  shaft  center  line. 

At  y  =90  degrees,  there  is  no  component  of  radial  flow, and  the 
loading  obtained  is  determined  primarily  by  the  twist  distribution  in  the 
blade.  Thus,  the  inboard  sections  have  a  positive  loading  and  the  out¬ 
board  sections  have  a  negative  loading.  On  looking  very  closely,  one 
can  see  a  small  discontinuity  in  loading  at  the  2.  5-  and  97.  5-percent  span 
station?  when  \p  is  90  degrees.  This  discontinuity  is  caused  by  the  tip 
of  the  blade  changing  from  a  trailing  to  a  leading  edge  (and  vice  versa 
for  the  root)  as  the  component  of  radial  flow  switches  from  a  positive 
to  a  negative  value.  As  (/  approaches  180  degrees,  the  loading  outboard  on 
the  blade  changes  from  negative  to  positive, with  the  positive  loadings  on 
the  inboard  sections  becoming  smallei.  At  ^  =  180  degrees,  a  load 
distribution  approaching  that  of  a  low-aspect-ratio  surface  is  obtained, 
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Figure  24.  RESPONSE  OF  STOPPED,  RIGID  ROTOR  TO  INPUT  OF  POSITIVE  SHAFT  TORQUE;  6C  =  0.138  RADIAN,  a,  =  0  RADIAN 
(a)  AERODYNAMIC  LOAD  AND  FORCE  TIME  HISTORIES. 


Figure  24.  RESPONSE  OF  STOPPED,  RIGID  ROTOR  TO  INPUT  OF  POSITIVE  SHAFT  TORQUE;  6C  =  0.138  RADIAN,  as-  0  RADIAN, 
(b)  BLADE  DISPLACEMENTS  AND  BENDING  MOMENTS;  SHAFT-VELOCITY  RESPONSE. 


• 

except  that  the  existing  value  of  $  -  5.  2  radians  per  second  causes 
positive  loadings  in  the  vicinity  of  midspan.  The  cause  of  the  high- 
frequency  change  in  loadings  occurring  between  /  =180  and  360 
degrees,  when  the  loadings  are  in  the  vicinity  of  zero,  is  not  clear. 

At  f  =  270  degrees,  however,  the  calculated  loadings  are  exactly 
what  one  would  expect  as  a  result  of  the  angular  velocity  of  the  rotor 
building  up  to  6.4  radians  per  second  at  this  point. 

Note  that  rotor  drag  is  small  when  the  rotor  is  in  the  vicinity 
of  ^  -  0  and  180  degrees  and  is  at  the  expected  peak  wher  \p  -  90  and 
270  degrees.  Note,  further,  that  shifts  in  the  chordwise  center  of  pres¬ 
sure  explain,  in  part,  the  torsional  response  (however  small)  which  is 
noted  in  Figure  24.  Starting  at  time  equal  to  zero,  the  initial  torsional 
response  is  due  to  fluctuations  in  loadings  on  the  leading  portion  of  the 
downstream  blade  with  the  c.  p.  of  these  loadings  being  around  midchord. 
The  larger  torsional  response  ensuing  at  t  -  0.  2  second  is  the  result  of 
the  sharp  change  in  loading  at  the  87.  5-percent  span  location.  (Here,  the 
contribution  to  the  generalized  force  in  torsion  is  much  greater  than  is 
produced  by  moments  at  inboard  stations.  )  Later,  at  f  in  the  vicinity 
of  210  degrees,  the  chordwise  c.  p.  is  moving  from  the  midchord  to  the 
three-quarter  chord,  and  the  small  fluctuations  in  loading  at  the  blade 
tip  provide  sufficient  moment  excitation  to  excite  the  torsional  mode 
of  the  blade  again. 
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CONCLUSIONS 


The  following  conclusions  are  drawn  as  a  result  of  this  study: 

1.  It  is  feasible  to  simulate  the  dynamics  of  a  rotor -wake 
system  by  means  of  numerical  integration  of  the 
applicable  equations  of  motion.  However,  it  would 

be  preferable  to  employ  a  hybrid  of  digital-  and 
analog -computing  means  to  implement  the  desired 
calculations  in  a  more  efficient  manner. 

2.  The  efforts  made  in  this  study  to  develop  an  aero¬ 
dynamic  model  to  generate  valid  aerodynamic  loadings, 
irrespective  of  whether  the  flow  is  primarily  tangential 
or  radial  with  respect  to  the  lifting  surface,  are  con¬ 
sidered  successful.  Additional  refinements  are  indi¬ 
cated  for  the  wake  model  in  order  to  represent  properly 
the  roll-up  of  wake  vorticity  trailed  under  conditions 

of  predominantly  radial  flow. 

3.  It  appears  that  no  aerodynamic  and/or  structural 
response  phenomena  of  unusual  nature  occur  during 
power  loss  and  braking  transients  for  the  rotors 
treated  in  this  study.  This  conclusion  could  pre¬ 
sumably  be  modified  by  performing  a  more  exten¬ 
sive  examination  of  the  bending  moment  distri¬ 
butions  along  the  span. 

4.  The  aerodynamic  loadings  and  structural  displace¬ 
ments  produced  on  a  rotor  stopped  at  0  -  180  degrees 
and  then  turned  through  one  rotor  revolution  are 
small  and  bounded. 


RECOMMENDATIONS 


1.  The  refinement  of  the  developed  digital -computing  program 
to  obtain  a  more  efficient  calculation  procedure  should  be 
performed  by  personnel  bearing  helicopter-design 
responsibilities. 

2.  There  is  need  for  extending  the  current  understanding  of 
the  aerodynamic  behavior  of  lifting  surfaces  having 
geometry  and  flow  orientation  typical  of  rotary  wings 

in  forward  flight.  This  extension  should  be  concerned 
with  both  steady  and  unsteady  flow  processes. 

3.  Calculations  should  be  made  to  determine  the  extent  tv. 
which  the  fulfillment  of  a  three-dimensional  boundary 
condition  and  the  employment  of  a  wake  model  would 
explain  discrepancies  existing  between  the  performance 
of  rotary  wings,  as  measured  and  as  predicted  by  calcu¬ 
lation  procedures  employing  "strip"  the  --ri  the  tra¬ 
ditional  assumption  of  uniform  inflow. 

4.  The  analysis  and  digital-computing  program  developed 
in  this  study  should  be  employed  to  investigate  problems 
of  rotor  stability  at  the  high  advance  ratios  that  will 
result  from  an  attempt  to  achieve  higher  forward  flight 
speeds.  Conclusions  drawn  by  others  relative  to 
instability  phenomena,  deriving  either  from  mechanical 
or  aerodynamic  sources  or  a  combination  of  these,  should 
be  put  to  the  test  with  the  assistance  of  the  developed 
program. 
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5. 


Efforts  should  be  initiated  to  augment  the  transient-blade- 
loads  program  for  purposes  of  examining  the  spacial  and 
temporal  distribution  of  loadings  and  stresses  encountered 
by  rotary  wings  upon  entry  into  a  vertical  gust. 
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APPENDIX  I 


THE  KINEMATICS  OF  ROTOR  BLADES  FREE  TO  FLAP,  LAG,  PITCH, 
BEND,  TWIST,  AND  ROTATE  WITH  A  VARIABLE  ANGULAR  VELOCITY 


The  objective  of  this  appendix  is  to  document  the  significant 
kinematic  relationships  that  are  utilized  in  (1)  performing  wake  calcu¬ 
lations,  (2)  solving  for  the  unknown  bound  circulations,  and  (3)  deriving 
the  equations  of  motion  that  yield  the  accelerations  in  each  assumed 
generalized  coordinate.  In  the  third  task,  the  evaluation  of  kinetic  and 
potential  energies  of  the  rotor  system  utilized  an  approach  that  was  based 
on  assumptions  (e.  g.  ,  the  location  of  the  elastic  axis,  the  axis  of  built-in 
twist,  etc.  )  that  were  slightly  different  from  those  adopted  to  compute  the 
displacement  and  velocities  of  points  relative  to  the  axis  system  de¬ 

fined  in  Figure  7.  These  different  assumptions  stemmed,  in  part,  from 
the  uncertainties  that  existed  with  respect  to  actual  construction  practice 
and,  in  part,  from  the  different  procedures  that  were  used  to  derive  kin¬ 
ematic  relationships  in  these  separate  instances. 

In  order  to  derive  the  *-y-f  position  (see  Figure  7  for  the  defi¬ 
nition  of  this  axis  system)  of  points  on  a  blade  as  defined  by  their 
location  (see  Figure  9),  use  was  made  of  the  axis -transformation  proce¬ 
dure  referred  to  in  Reference  6  and  outlined  in  more  detail  in  Reference  21. 
The  procedure  is  implemented  by  starting  with  a  point  on  the  blade  and 
moving  inboard  to  the  hub  by  successive  translations  and  rotations  based 
on  all  of  the  generalized  displacements  assumed  to  exist  in  deriving  the 
equations  of  motion  of  the  rotor.  The  degrees  of  freedom  assumed  for 
this  analysis  require  the  following  transformations: 

1.  Rotation  of  the  blade  about  the  elastic  axis  to 
eliminate  the  twist  due  to  torsion. 
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2.  Translation  of  the  blade  to  remove  bending  displace¬ 
ments  in  the  flapwise  and  chordwise  directions 
(shortening  of  the  blade  in  /  direction  is  neglected). 

3.  Translation  of  the  elastic  axis  to  the  feathering  axis. 

4.  Rotation  of  the  blade  about  the  feathering  axis  to 

>!< 

eliminate  the  built-in  twist. 

5.  Translation  of  the  ^  axis  to  that  position  along  the 
blade  span  assumed  to  be  the  origin  of  the  £  axis. 

6.  Rotation  of  the  root  chord  about  the  feathering  axis 
to  a  plane  parallel  to  the  shaft  plane. 

7.  Translation  of  the  origin  of  the  feathering  axis  to  the 
drag  hinge. 

8.  Rotation  about  the  drag  hinge  such  that  the  feathering 
axis  is  parallel  to  the  radial  direction. 

9.  Translation  of  the  origin  to  the  flap  hinge. 

10.  Rotation  about  the  flap  hinge  such  that  the  feathering 
axis  is  parallel  to  the  shaft  plane. 

11.  Translation  of  the  origin  from  the  flap  hinge  to  the 
center  of  the  rotor  shaft. 


It  should  be  recalled  that  this  assumption  is  incorrect.  The  practical 
consequences,  however,  are  not  believed  to  be  significant. 
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12.  Rotation  of  the  axes  (established  by  the  previous 

transformation)  about  the  rotor  shaft  to  line  up  with 
the  nonrotating  axis  system  defined  in 

Figure  7. 

13.  Rotation  of  the  x,  -  y,  ■  y.  axes  about  the  y,axis  to 
obtain  fhe  "wind-axis"  system. 

On  making  the  usual  small-angle  assumptions,  and  on  dropping 
second -order  terms,  the  following  general  expressions  for  blade  coor¬ 
dinates  are  yielded  by  the  above -outlined  transformation  procedure; 
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These  expressions  are  utilized  to  determine  the  coordinates  of 
control  points  and  the  end  points  of  vortex  elements,  both  fixed  and  free, 
in  the  x-y-}  space.  For  example,  the  x-y-f  location  o f  the  three -quarter 
chord  point  located  in  blade  t  at  the  center  of  blade  segment  j  and  at 
time  m  is  given  by 
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where 


A 

*.  J 

6m,  i  *  6j 

rj-  -  q  distance  from  elastic  axis  to  three-quarter 
chord  at  center  of  blade  segment  J  . 

-  lumped  value  of  6^,  i.  e.  ,  value  of  built-in 
twist  at  center  of  segment  J  relative  to 
chord  at  f  =  0. 

-  blade  twist  due  to  torsion  at  time  m,  blade  i , 
station  j  . 

Qm  L  -  pitch  angle  of  blade  chord,  at  £  =  0,  with 
respect  to  the  shaft  plane;  blade  i  at 
time  . 

The  coordinates  of  the  wake  are  defined  in  a  fashion  similar  to  that  out¬ 
lined  in  the  appendix  of  Reference  [12].  Since  the  rotor  in  TBL-II  rotates 
with  a  variable  angular  velocity,  the  time  increment  used  for  convecting 
the  wake  is  t--tm  .rather  than  ,  as  was  used  in  TBL-I. 

m  m.f  j} 

In  addition  to  determining  the  coordinates  of  the  blade  in  *-y-  y 
space,  it  is  necessary  to  determine  the  velocities  of  given  points  on  the 
blade  with  respect  to  inertial  space.  Since  these  velocities  are  used  to 
compute  the  angles  of  attack  used  in  the  "strip"  analysis  and  to  obtain 
the  velocity  component  normal  to  the  blade  when  satisfying  the  "three- 
dimensional"  boundary  condition,  it  is  more  convenient  to  work  in  the 
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shaft-axis  system.  Coordinates  yi(i,7)>  anc*  ^(£.7}  are 

thus  differentiated  to  obtain  )»  and  /p). 

If  the  radial,  tangential,  and  normal  velocities  of  a  point  on  a 
blade  relative  to  still  air  are  designated  as  Uj,  ,  and  r{t  and  C 

are  axes  located  in  a  plane  tangent  to  the  blade  surface  at  the  selected 
point),  these  velocities  are  given  by 

U'f 

(10) 
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where  e  %  and  are  unit  vectors  oriented  along  the  respective  f,  rj, 
and  f  axes,  and  ^  is  the  total  velocity  vector  of  a  point  on  the  blade  as 
defined  in  the  shaft-axis  system,  lfu.tr,  and  i?are  designated 

as  components  of  the  velocity  that  is  induced  by  the  rotor  wake  in  z-y-y 
space  and  as  is  taken  as  a  small  angle,  Figure  7  shows  that 
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where  V  ia  the  forward  velocity  of  the  hub  along  the  negative  x  axis.  Equa¬ 
tion  (10)  can  be  employed  to  obtain  the  desired  velocity  components  by 
noting  that 
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where 


On  differentiating  the  coordinates  t//(4.ff)  and  both  with 

respect  to  time  and  with  respect  to  §  and  y,  and  on  neglecting  second- 
order  terms,  the  following  expressions  are  obtained  for  (J^ ,  6^  ,  and  , 
where  these  unprimed  components  designate  air  velocities  relative  to  the 
blade  rather  than  blade  velocities  relative  to  still  air: 
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Note  that  the  above  velocity  components  have  been  expressed  in 
terms  of  an  ^-^•{'set  of  axes  whose  origin  is  located  at  an  i\j  control  point. 
If  rj.  is  assumed  to  be  a  distance  from  the  elastic  axis  to  a  three-quarter 
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chord  point,  then,  for  the  "strip"  analysis,  the  angle  of  attack  at  the  three- 
quarter  chord  can  be  expressed  as 


£/„  . 

or  •  =  tan'  t*L±dL  . 

m>  *>  j  (j 

The  boundary  condition  for  the  "three-dimensional"  analysis  is  giver  by 


where 


j 


0 


0m0 

-  velocity  induced  by  the  bound  circulation  on 
blade  i  at  the  i,j  control  point  . 


The  elastic -blade  degrees  of  freedom,  v'ttu',  and  p,  are  assumed 
to  be  represented  by  a  superposition  of  assumed-mode  shapes.  Combined 
flapwise-  and  chordwise -bending  modes  can  be  represented  as 


where 


<?>  ,(»  M  .1 


^Jare  the  flapwise-  and  chordwise-displacement  mode 
shapes,  respectively,  for  the  p**  mode. 

k(*  is  the  ratio  of  the  chordwise  displacement  of  the  blade  tip 
to  the  flapwise  displacement  of  the  blade  tip  for  the  p*h  mode 


and 


is  the  flapwise  displacement  of  the  blade  tip  in  the  p*h  mode. 


Blade  torsion  is  represented  as 


APPENDIX  II 


AN  AERODYNAMIC  MODEL  FOR  CALCULATING  LOADINGS  ON 
ROTOR  BLADES  OPERATING  AT  HIGH  ADVANCE  RATIOS 


The  horseshoe-vortex  representation  employed  in  TBL-I  con¬ 
stituted,  in  essence,  a  classical  lifting-line  approach  tor  evaluating  the 
spanwise  distribution  of  circulation  on  a  high-aspect-ratio  surface.  In 
accordance  with  the  Prandtl  simplification,  the  lifting  line  located  in  a 
given  spanwise  segment  of  the  blade  is  considered  to  extend  to  plus  and 
minus  infinity  in  stating  the  boundary  condition  to  be  satisfied  at  the  three- 
quarter  chord  of  that  same  segment.  The  net  result  of  this  approximation 
is  that  the  loading  on  a  given  spanwise  segment  is  given  by  the  lift  and  drag 
of  a  two-dimensional  airfoil  whose  angle  of  attack  is  specified  by  the  flow 
velocities  required  to  achieve  flow  tangency  at  the  control  point  of  the 
segment  under  consideration.  As  was  discussed  in  the  report,  it  follows 
that  the  applicable  angle  of  attack  of  the  two-dimensional  segment  is  given 
by  the  kinematic  properties  of  the  blade  and  the  downwash  induced  both  by 
the  free  wake  and  the  bound  circulation  in  blades  other  than  the  blade  in 
which  the  segment  is  located. 

The  downwash  caused  by  elements  of  bound  vorticity  located  in 
blades  other  than  the  blade  being  considered  was  neglected  in  the  calcu¬ 
lations  performed  in  TBL-I.  A  subsequent  introduction  of  these  elements 
of  downwash  into  the  computational  procedure  demonstrated  that  these 
elements  are  a  negligible  contribution  at  low  advance  ration.  Even  though 
these  contributions  may  remain  small  at  high  idvance  ratios,  these  effects 
have  been  retained  in  the  TBL-II  model,  in  view  of  the  requirement  for 
computing  the  velocities  caused  by  the  elements  of  bound  vorticity  located 
in  the  blade  containing  the  control  point  at  which  the  boundary  condition 
is  being  satisfied.  This  latter  requirement  stems,  of  course,  from  the 
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elimination  of  the  Prandtl  simplification  (to  be  discussed  below),  since 
this  simplification  is  permissible  only  when  calculating  loadings  on 
large-aspect-ratio  surfaces. 

As  was  outlined  in  the  body  of  the  report,  the  objective  of  the 
TBL-II  study  requires  the  computation  of  loadings  on  lifting  surfaces 
exposed  to  flows  of  varying  obliquity.  The  1.  licopter  rotor  blade  at 
a  high  advance  ratio  encounters  flow  obliquities  which  cause  it  to  change 
continuously  from  a  high -aspect-ratio  surface  to  a  low-aspect-ratio 
surface  as  it  changes  in  azimuth  through  90  degrees.  However,  it 
was  recognized  that  the  development  of  an  aerodynamic -loading  theory 
for  a  surface  exposed  to  the  flow  alignments  encountered  by  a  slowly 
rotating  rotary  wing  is  a  research  task  of  considerable  magnitude. 
Accordingly,  a  limited  objective  was  selected,  namely,  to  develop 
a  loading  ■' calculation  procedure  utilizing  as  simple  a  vortex  model 
a.3  possible,  where  the  validity  of  this  model  would  be  judged  by  its 
ability  to  yield  correct  loadings  at  the  two  extremes  of  low  and  in¬ 
finite  advance  ratio.  As  discussed,  in  part,  in  the  report,  the  achieve¬ 
ment  of  this  objective  called  for  the  following  modifications  to  the 
aerodynamic  model  described  in  Reference  [12]: 

1.  Move  the  spanwise -oriented  bound  vortex  from  the 
one-quarter  chord  to  the  leading  edge  in  order  to 
obtain  a  suitable  vortex  geometry  for  flows  that  are 
highly  radial  with  respect  to  the  blade. 

2.  Place  the  control  point  at  the  center  (50-percent 
chord)  of  the  blade  segment. 

3.  Introduce  tests  to  determine  whether  root  or  tip 
of  blade  is  a  leading  or  trailing  edge. 

4.  Introduce  logic  that  will  trail  vorticity  from  the 
rounded  "leading"  edge  of  the  blade  in  the  reversed 
flow  region. 
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5.  Add  effects  of  spanwise  and  chordwise  bound 
vorticity  (from  all  blades)  to  the  matrix  of 
induced-velocity  coefficients. 

6.  Redefine  equations  yielding  the  unknown  bound 
circulations  in  terms  of  the  flow-tangency  re¬ 
quirement  at  the  selected  control  point. 

7.  Develop  methods  for  (a)  adding  viscous  forces 
to  the  circulatory  forces,  (b)  approximating  the 
effects  of  Mach  number,  and  (c)  estimating  the 
tangential  location  of  the  center  of  pressure. 

Modifications  (1)  through  (4)  rre  straightforward  in  concept 
and  execution.  They  involve  changes  in  kinematic  definitions  having 
an  influence  on  wake  geometry  and  the  motion  of  the  control  point. 
Accordingly,  these  changes  are  not  documented  in  this  appendix. 

Modifications  (5)  and  (6)  are  key  items  which  convert  the 
load-prediction  scheme  lor  high -aspect-ratio  surfaces  to  one  that 
will  predict  loadings  also  for  low-aspect-ratio  surfaces,  as  per 
the  example  case  discussed  in  the  body  of  the  report.  The  method 
employed  in  TBL.-II  to  introduce  the  effects  of  bound  vorticity  into 
the  downwasb  computation  is  ^ive.i  below. 

If  an  induced-velocity  coefficient  is  defined  to  represent  the 
combined  influence  of  a  hoi seshoe -vortex  element  fixed  in  the  blade  and 
the  adjacent  elements  of  free  wake  (assumed  to  have  the  same  strength 
as  the  bound  vortex),  the  resulting  composite  velocity  coefficient  can 
be  written  as 


where 


°~ik  2  1  ~*rik‘lA  +  1  are  induced -velocity  coefficients 

of  the  adjacent  free  wake. 

_  m,  i,  J 

<7^  j  is  the  induced-velocity  coefficient  of  a  spanwise- 
oriented  element  of  bound  vorticity. 

are  induced  velocity  coefficients  of  the 
chordwise-oriented  elements  of  bound  vorticity. 

In  the  "strip"  analysis  applicable  to  low  advance  ratios,  the 
Prandtl  simplification  is  implemented  by  (a)  neglecting  all  contri¬ 
butions  from  the  barred  coefficients  when  ib  =  i  ,  and  (b)  assuming 

—  /"i  A 

that  us.  #  =•  _ — £-  .  In  this  manner,  it  becomes  possible  to  compute 
V  2lrbs 

section  angles  of  attack  and,  thus,  to  use  experimental  lift  and  drag 
data.  In  a  three-dimensional  analysis,  the  velocity  induced  at  a  con¬ 
trol  point  by  the  ib  *  i  components  of  bound  vorticity  approaches  that 
induced  by  a  constant  strength,  doubly  infinite  lifting  line  as  the  flow 
becomes  primarily  tangential  at  low  advance  ratios.  Thus,  at  low 
advance  ratios,  a  loading  result  roughly  equivalent  to  that  yielded 
by  the  "strip"  analysis  would  be  expected.  (In  effect,  the  retention 
of  the  tb  z  c  components  of  bound  vorticity  accounts  for  the  influence 
of  finite  wing  length. ) 

Notwithstanding  this  qualitative  agreement  between  the  "strip" 
and  the  more  exact  analysis  at  low  advance  ratios,  significant  differ¬ 
ences  would  be  obtained  in  actual  calculations  since  one  method  uses 
the  lift-curve  slope  of  an  actual  airfoil  immersed  in  a  fluid  of  finite 
viscosity  and  the  second  method  uses  the  theoretical  iift-curve  slope 
of  a  flat  plate.  There  would  also  be  differences  in  calculated  loadings 
as  a  result  of  variations  in  experimental  lift-curve  slope  with  Mach 
number. 
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The  calculation  procedure  developed  for  low-aspect-ratio 
surfaces  was  arbitrarily  modified  such  that,  at  low  advance  ratios, 
calculated  loadings  approximate  those  yielded  by  the  "strip"  analysis. 
This  modification  consists  of  altering  the  strength  of  the  induced-velocity 
coefficient  of  the  tb  -  £  bound  vortex, as  a  means  of  representing  the 
circulation  existing  on  a  real  airfoil, in  contrast  to  the  circulation  that 
represents  the  bound-vortex  equivalent  of  the  circulation  distributed 
across  the  surface  of  a  flat  plate.  The  appropriate  correction  to 
induced-velocity  coefficients  associated  with  radially  or  tangentially 
oriented  components  of  bound  vorticity  is,  for  example,  given  by 

_  m,  i $  j  2ir  _  m,t,j 

^ * *. *  J dcL\  C7T^-«./ 

Jmj 

dLCc  \ 

Note  that  ^mi  , m  can  be  obtained  from  experiment  as  a  function 

of  Mach  number,  a  value  of  which  can  be  computed  for  every  <  j  con¬ 
trol  point. 

The  appendix  in  Reference  [12]  shows  that  the  unknown  bound 
circulations  required  to  satisfy  the  boundary  condition  (as  expressed 
in  the  "strip"  analysis)  satisfy  the  following  equation  set: 

</] 


The  column  matrix  ym  (  j  |  represents  the  unknown  circulation  strengths, 
and  the  column  matrix  {  r*  l  represents  that  portion  of  the  circulation 
at  each  spanwise  segment  as  determined  by  all  ge.  metric  and  downwash 
contributions  to  the  local  angle  of  attack  other  than  the  downwash  created 
by  the  adjacent  elements  of  free  wake  and  the  elements  of  bound  vorticii.y 


in  all  blades  other  than  blade  cb*i.  The  square  matrix,  c-l-J 

L  ‘  A.  /_ 


,  represents 
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the  composite  influence  of  the  induced-velocity  coefficients  of  the  free-wake 
elements  adjacent  to  all  blades  and  the  bound  vorticity  in  all  blades  other 
than  the  blade  containing  the  i,j  control  point.  Each  term  in  the  square 
matrix  is  given  by 


°7  6.i 


-b: 

* 


cLcc/m.  i,  j 


_ J  ,  — 

U4i 


+  <rib.t  ~ 

Ib/i  U  4  i 


(11) 


except  for  the  terms  on  the  diagonal  which  are  given  by 

dCL 


°£U.j 


da 


j  *- 


__  m-*,4 

~ib  •  i,  2. 


<n*, j 


J *  • 


The  implementation  of  item  (6)  (see  above)  requires  that  the 
err  i. ion  set  as  derived  for  the  "strip"  analysis  be  reformulated  in  terms 
nree-dimensional  boundary  condition.  Note  that  the  condition  of 
flow  tangency  at  a  control  point  may  be  expressed  as 


.  Q  +*+ 

UJ  +  cos  &m 


(12) 


where  i/1  is  the  velocity  of  air  normal  to  an  it\  point  on  the  blade 

as  derives  from  geometric  and  kinematic  considerations  only,  and  olt^  ^  ■ 

is  the  total  ^-component  of  downwash  created  at  point  t,  j  by  all  fixed  and 

free  components  of  circulation  in  the  rotor -wake  system.  Since  the  total 

downwash,  <2^,  »  consists  of  the  sum  of  the  downwash  created  by  (a)  all 

elements  of  the  wake  except  for  the  adjacent  wake,  urT  ,  and  (b)  the 

adjacent  wake  plus  all  bound  vorticity  in  the  blades,  u ^  ,  Equation 

*,  J 

(l<d)  can  be  expressed  as 


^4 


‘r 

cos  &L 


- 


r  •  . 


(13) 
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The  downwash  component,  u^A  ,  may  be  further  subdivided 

*,  J 

into  (a)  downwash  caused  by  the  adjacent  wake  elements  and  the  bound 
vorticity  on  blades  other  than  that  containing  the  i,j  control  point,  and 
(b)  downwash  caused  by  bound  vorticity  on  the  blade  containing  the  t,J 

control  point.  Accordingly,  we  find  that  the  quantity  urA  ,  may  be 

m,  it  J 

expressed  as 


/V£  iet/ 


t*  r-  r  ,j  i  LQ  r 

J  ^  * 

The  O“0\‘J  matrix  in  the  above  expression  is  identical  to  the  quantity 

*-  .4  r  •  -| 

found  within  the  brackets  in  Equation  (11).  The  ; .  .  matrix  con- 

L  »*  ‘i.i  J 

tains  that  portion  of  the  original  definition  of  ^Tq"1  ,  which  is  omitted  in 

*k4 

the  formulation  of  the  "strip"  analysis  i  tus,  the  three-dimensional 
boundary  condition  may  be  expressed  as 


where 

This  form  of  the  boundary  condition  shows  that  it  is  only  necessary  to 
evaluate  some  additional  (^coefficients  in  order  to  modify  the  "strip" 
analysis,  as  per  the  objective  stated  in  Item  (6), 


As  was  indicated  in  the  body  of  the  report,  a  "two-dimensional" 
stall  limit  is  imposed  on  the  values  of  circulation  that  are  obtained  from 
a  solution  of  Equation  (14).  This  limiting  procedure  restricts  the  strength 
of  the  tip  vortex  in  the  wake  to  values  compatible  with  the  maximum  circu¬ 
lation  that  can  be  produced  by  a  stalled  blade.  Thus,  the  boundary  condition 
is  not  satisfied  at  all  control  points  if  and  when  a  portion  of  the  blade  is 
considered  to  be  stalled  at  low  advance  ratios. 
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Application  of  the  Kutta - Joukowaki  law  to  obtain  aerodynamic 
loadings  from  the  computed  bound  circulations  yields  forces  on  the  blade 
that  derive  only  from  potential  flow  considerations.  Item  (7)  was  thus 
implemented  by  arbitrarily  introducing  a  tangential  viscous  force  based 
on  the  two-dimensional  angle  of  attack  that  is  approached  at  low  advance 
ratios.  When  the  flow,  relative  to  the  blade,  is  primarily  tangential  in 
direction,  the  normal  force  on  the  blade  is  given  by 


L5_ 


>.*.  j 


and  the  tangential  force  (chordwise)  is  given  by  a  circulatory  component. 


DS  *  -/>  Ur 
and  a  viscous  component, 


°  ^  )  2  bj  C°7 m. i 


to 


where  the  drag  coefficient,  CD  ,  is  determined  as  a  function  of  as 

defined  in  the  body  of  the  report.  Radial  components  of  flow,  .  .  ,  also 
produce  a  normal  force  on  the  blade  which,  because  of  the  lumping  process 
used,  are  lumped  forces  at  the  radial  locations  of  tangentially  directed 
elements  of  bound  vorticity.  These  concentrated  loads  are  given  by 

Ufm  t  j\fm.  ‘  ] 

which  loads  are  then  "smeared"  to  yield  a  radial  distribution  of  loading. 


Pitching  moments  about  the  feathering  and  elastic  axes  are 
assumed  to  derive  from  quasi-static  lift  acting  at  the  quarter  chord  and 
from  the  quasi-static  moment  about  the  quarter  chord  resulting  from  angular 
velocities  in  pitch  and  twist.  The  center  of  pressure  of  the  quasi-static 
lift  caused  by  pure  effective  plunging  is  also  assumed  to  shift  rearward 
as  a  function  of  the  local  two-dimensional  angle  of  attack  and  the  local 
Mach  number.  Curve  fitting  of  experimental  data  obtained  with  NACA 
0012  and  0015  airfoil  sections  yielded  the  following  expressions  for  the 
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APPENDIX  III 


THE  DERIVATION  OF  THE  EQUATIONS  OF 
MOTION  NEEDED  IN  THE  TBL-II  STUDY 


A  Lagrangian  approach  was  used  to  derive  the  equations  of  motion 
for  a  rotor  with  articulated,  flexible  blades.  As  was  discussed  in  the  body 
of  the  report,  the  introduction  of  the  appropriate  kinematic  constraints 
and  the  selection  of  the  proper  antisymmetric  and/or  symmetric  bending 
modes  convert  equations  applicable  to  an  articulated  rotor  to  equations 
applicable  to  a  teetering  rotor  or  the  so-called  "rigid"  rotor. 

As  a  result  of  specifying  (1)  the  existence  of  offsets  (as  defined 
in  Figure  8)  and  (2)  the  blade  to  flap,  lead-lag,  and  feather  in  the  order 
given,  the  following  expressions  are  obtained  for  the  x.-y-j  components 
(see  Figure  9)  of  the  total  absolute  velocity  of  a  mass  point,  P,  located 
in  a  blade  by  the  rj  coordinates  defined  in  Figure  9.  (It  should  be  noted 
that  the  blade  cross  section  is  assumed  to  be  thin  and  symmetrical.  ) 

\  *  i  +  u  +0)  -  ry 
Vi  *  y  +  V  +  rx  -  /?  J. 

V}=  j  *  pi  -  fx 

where 

x  «  -i 

(/  -  if  COC  G^  -  ufssn  -  0ystn(df  +  4>) 

|  >  ur  cos  9 i  t  if  sin  9*  *  c-os  ( + 4) 

t  *  i  -S 

y  »  ea  *  if  cos  Gf  -  ur 

j  =  ur  cos  Qj  +  %r sin  0f  *■  y  sin(0f  *  <f>) 

u  *  V  {(Cp  -  c„  sm/S)  sin  C  *(eD  cos  +  e„  sin  c  -  ip)  cos/S]  -jiCp  cos  c~  Cep 

V  *  <i>  {[ tfcos{  +  (e0  co*(-c0  s in{  +  ep)  cos  fi  -  Zp  sin 3  cosf\  cos  & 
cos  {  -  epSin c]swfi sin  &}  -  fep  cos  G 
•fl  I Cp  sin  C  CoS  6  -  (eD  *  ep  COS  f  -  ep  sin  £)  sin  0  j 
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*  '  *  {l*  d~(%  cos  (  t-e^  sin  C)]  sin  ft  cos  6 

-  \tfi  cos  C  +(e0  cos  C  -  «#  sin  C  *■  C?)  cos  /S  -Cp  sinfi  cos  stn&J 
+-/4{fes  +  ef,  cos  C  -  Cp  sinf)  cos  6  +cf  sin  £  sin  0  }  «■  C  &p  Jin  9 

p  *  <f>  Sen cos  C  +  Jin  C  & 

y  »  <js  sen  fi  sen  fi  cos  6  +  cos  fi  stn  9  J  -fi  cos  fi  cos  6  -  fi  sin  & 
r  *  f  f  cos  fi  cos  9  -sen  fi  sin  fi sin  6  j  + fi  Cos  fi sen  6  -  Ceos  0 

and  <5  is  the  shortening  in  the  ^-direction  due  to  the  blade  being  bent. 

The  shortening,  $  ,  oi  the  blade  may  be  expressed  as 
S'A  *■  ysert(&f  +  0)  +  ^  rjc.o*(  e€  +-  4>) 

where 

or,  »  ur  cos  +  irsin  9^ 
l/'f  »  iA  COS  Oj.  *  isf'sin  9^ 

and  H  is  the  shortening  of  the  elastic  axis  (in  the  *-dir ection)  up  to 
station  f  and  is  given,  approximately,  by 
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In  order  to  express  the  kinetic  and  potential  energy  of  a  single 
blade,  it  is  convenient  to  expand  the  blade  displacements  i rt  u<tnd  f>  in 
terms  of  the  generalized  coordinates,  and  <p(V ,  defined  in  Appendix  I. 

The  kinetic  energy  of  a  single  blade  is  given  by 


The  potential  energy  can  be  expressed,  following  Reference  [19],  as  a 
function  of  the  incremental  strain  energy  caused  by  elastic  bending  and 
torsional  deformations  of  a  thin  twisted  blade  under  initial  tension.  This 
initial  tension  is  a  function,  primarily,  of  the  shaft  rotation  velocity,  tjf  . 
Thus,  the  incremental  strain  energy,  1C  ,  is  given  by 


where 


U  •  {j'uLi,  (‘*£1,  V) if  <■  \l\oJ^e,(TPi 
-/*XeAXdf 


initial  tension  on  blade  section  at  station  (  and,  on 
ignoring  contributions  due  to  lagging  and  flapping, 
may  be  expressed,  approximately,  as: 

.  r* 


aL 


i>  J  («n  *  f) 

dd.  die.  fdOz\z 

of  Tf  + 


tr,  -- 

fr,  = 


= 


dlir  _  ,  d9t  cLur  <**6,  fdOA1. 

dp  2  ^  3T  sd  "  '  U77  ’ 


d£ *  df  d{  d(r-  \~zii 

flexural  rigidity  about  major  neutral  axis 
flexural  rigidity  about  minor  neutral  axis 
torsional  rigidity 


*  =  ilidA 


ill 


li  =  t/  l‘dA 
B,  •  l  ftf-k',)** 

6*  ~  JA  7  ('?*’-  k*}  *■* 

and  A  is  the  blade  cross-section  area  effective  in  taking  tension,  with  y 
being  the  distance  from  the  elastic  axis  to  the  element,  dA  (positive 
when  dA  is  closer  to  the  leading  edge). 


On  applying  Lagrange's  equation  to  obtain  the  equations  of  motion 
for  a  rotor  translating  with  a  fixed  velocity,  with  the  generalized  coordinates 
considered  to  be  and <p(V  ,  one  finds  that  the  required  partial 

differentiation  and  the  integration  over  the  mass  of  the  blade  lead  to  a 
multiplicity  of  terms,  many  of  which  are  of  second  and  higher  order.  In 
this  study,  higher -order  terms  were  dropped  and  small-angle  assumptions 
were  made,  when  permissible. 


For  example,  we  show  below  the  equation  that  results  for  a  blade - 
bending  degree  of  freedom,  namely,  the  r  mode  of  coupled  flatwise  and 
chordwise  bending: 


O **£?+&* 

F  f 

F  i  J 


ir> 

•ur 


c?  us*' 


=  &.F. 


On  defining  the  sectional  properties  of  the  blade  as 
m  -  mass  per  unit  span,  slugs/ft 

<•  =  chordwise  distance  of  c.  g.  from  the  elastic  axis,  ft 

^  =  mass  moment  of  inertia  about  the  elastic  axis  per 

unit  span,  lb-sec^ 


and  on  noting  that  the  blade -bending  mode  shapes,  as  derived,  are  orthog¬ 
onal  (including  flapping  which  can  be  considered  to  be  a  zero-order  bending 
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node),  we  find  that  the  above  coefficients  are  given  by  the  following  integrals 


A<r) 

M  <Puf  * 

J  r»(cn  +f)  | 

!  -  *sin  (8t 

+  9)  +■  itv  Vj.  1  cos  (&t  +  9) 

\dU 

A<r  , 

H/Sus 

O 

4?>  , 
*fV 

f 

Jo 

m(epr£)  j 

(&(  + 

Z3  )  /  (r)  C  (r> 

0)  ~  k.v  +„.  cos 

(Of  +o)\ 

cL£ 

A<r> 
A0ur  ' 

J  m  |  e0  sin  Qf 

cc 

+  (*o 

cos  0(  +  e) 

\dt 

c\('r)  = 

LVOV' 

o, 

r  /  P 

A(r‘r)  - 

’’  cV  oV  ~~ 

/ 

f  [c 

4."'] 

> 

uo 

d'J?  =  O 


Jo  m  scn(ef  +  9)  cos  (O'  +&)d? 

~  f  nz( 6f  *■  O) aL£ 

-  f  ^  cos  z  (0f  +■  &)  d( 

V  o 

where  g 

A(()  *  [  m(c#+-f)d€ 

..  (r)  _  df£}  d&f  (r)  (r) 

and  *  <**  df  *  " 

,  (r)  ,(r)  dfr  d0.  (r> 

v  -  k  — —  ■& - L  / 

dt  dt  " 
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c'S?  =  f  *m  (e*  *  f) 

f  J  }Jin  *(6f  *  9)  d  £ 

+  f*r»e  /^klr>  W'str,(ef  *  6)  cos  (0g  +  0)cL{ 


#u 

dor 


(r) 


-z*-™  Fjr-1”  -z^cc?’  -  f. 


<r) 

US’ 


r 

where 


Cr,''-J'{ei,  V'V,"'*  et,  /f  V»"’] 


W-J,  « 


F~’m  r‘/'A(<)e4/i,r)‘*{ [fo*8U)cAtir)<d(+/g*A(t)eAfr>dt]. 
(The  integrals  yielding  the  generalized  force,  o  F<r)  ,  are  discussed  and 
defined  in  Appendix  IV. ) 


In  the  above  example,  all  terms  have  been  retained  except  those 
inertial  terms  that  become  zero  as  a  result  of  assuming  orthogonal  modes. 
Although  the  inertial  coupling  between  bending  and  angular  acceleration  of 
the  rotor  shaft  is  small,  these  terms  have  been  retained  as  a  result  of  f 
being  very  large  whenever  a  braking  torque  is  applied  to  the  shaft  of  the 
rotor.  A  similar  consideration  applies  to  the  retention  of  the  inertial 
coupling  between  flapping  and  shaft  acceleration. 


In  setting  up  the  equations  of  motion  in  matrix  form  for  machine 
solution,  it  is  advantageous  to  expand  further  certain  integrals  to  separate 
them  into  a  time-independent  and  time-dependent  part.  Tables  III  and 
IV  present,  respectively,  the  equations  of  motion  developed  for  this  study 
and  the  coefficients  required  to  repres,  .t  the  inertial,  elastic,  and  centrif¬ 
ugal  properties  of  a  rotating,  twisted  blade  experiencing  large  feathering 
displacements.  It  should  be  noted  that  Coriolis  acceleration  terms  have 
been  retained  only  in  the  equations  for  the  shaft-rotation,  flapping,  and 
lead-lag  degrees  of  freedom. 
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TABLE  III 

EQUATIONS  OF  MOTION  FOR  AN  ARTICULATED 
ROTOR  WITH  NB  FLEXIBLE  BLADES 


Shaft -Rotation  Degree  of  Freedom: 


N6'iA„  +£  *  2 f  [2  A  (5f,  3„)] 

i+f  l<i/  4  ' J 

Flapping  Degree  of  Freedom: 

+  A  *•>  *  44»,  '■21  A r  ^ 

+  *‘[b„+A  bm] 

Lead-Lag  Degree  of  Freedom: 

**  V  6  V  *  *<«.■ *  Z  *  z 

p  ‘  f 

Feathering  Degree  of  Freedom 

^  'A  V  y  <■><»,  y«i-4„  *■!"*}£■ 

y  rfc *BuA +Z  [f«y 

Coupled  Flapwise  and  Chordwise  Bending  Degree  of  Freedom 

v<~.  *$*&+$*%  *  *r*~+n? c" 

*  *  27^/"^;'’/  *  -.'V  C;-Z*?'W  - «  y 

C  *  ‘  ^  <  J  *  j 

Torsion  Degree  of  Freedom 

VAP0i  +A  A**i  ’h?iAf4.  *  &i  A0t  +  Z.  ^  ^AA00 

P 

+  (s>f)  A.  Y /A(*>rii'f)^  r<f'l  <'>s-(P-*\  cft}  , 

^  (Z.^;  ^  j  -2.^t  ^  ^f0.  :0F 
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TABLE  H 

DEFINITION  OF  COEFFICIENTS  (INTEGRALS) 
APPEARING  IN  TABLE  in 


r*  t 
Aff  -  J  m(c.n  +f)  cL( 

An  *  ~  J  rr>(c*  +f)(eF  *■  {)<*■{ 

Bpt.  -  m(t#i-4)Cp  cL$  +  sin  &t J  m  (e*  *■  f)(e0  +  e)d$ 

+  cos0i  J'  m(cHi-Z)z6{  d€\ 

o  J 

Bfi  -  -  f  mfc*  +  i)(t,  +2) 

uo 

&rr.  =  -  eof  (er  *4)  m(£{  -  ef  (tF  +  eD)  J  md£ 

O  0 

+  COS  0t(e,+eD) J  m(c.+  e)eL£  -sin  y  eD) J  me  9(d< 

fi,¥-  -(tf  +  eD) f  rr>(ep+<)d4 

Ar*t  ■  A'  \(*o '%)  J m(et^4)d{ +cosd;J m (e0+  c)(e,  * 4)d{  -sind^f mc6i(et+4)d(^ 

r*  ~  f*  f* 

'  Ai ‘•if  m(e>  r4)(ep*4'd{  +  cp(e0-er)  j  md(  +  cos6tZpJ  m(ee+e)d4 

Jo  Jo  V 

i  [cos 6^(t0-%) -sin 9- cp J /  me O^df  &  sin 9i (e0 -ep) j  m(e0  *  *)df 

+  cosZ6t;[  meet0.d4  *  sm8,  cos&- f  me^(e0  *■£<:)  d4  -bcos  26t  f /u9,  d( 

*4  *  Jo  Jo  f 

*  COS  sin  e*f0  ~coi6‘  stn  6if  {cpJm(cp+f)d4 

t- sin  9l J me.(c,+(Ut  •bcosQ^J  medf(e.p  +  t)di+ sin6;Jme(ep+Z)d.z\> 


+  cos 
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TABLE  JX  -  Continued 


fi/s 


1  /  m(e,+()*dt 

Jo 


A00(  m  cos6tJ*m(f',tt)(c0  *€)ctf  -sirt&ij  m(e,  *  f)  e0f  dt 
A/J’*  cos 0t  J*m(e, *() tf -sin6cJ*m(e,  +f)ef?)6sdt 


Bfi>i  ‘  ~Bf'i 


ffj$t  •  ~  3ft 

Cfil.  -  J m(ep*i)d(+sin6i^ ">(«,+ f)(e,+*)d{  +  COS 0-  J  m(cp+f)e&f 

Cgt  •  -f  m  (ep +{)(*,  +  ()  ct{ 

Jo 

Af£  -  ~ [  m(c„  *-Jp  (tp  +  {)d{ 

4  "o 

f* 

Acc  -  J  m(t,  +  f)*d( 

*  / * 
m(tg  +f)e  &fd{  +  Sin  0; J  rr>(ep  +  f)(e0  +  e)cL( 

A(%=  cos0^*m(e0u)^efdf  - / 

+Stn  &4-  [//">«,  +{)CJctt+f'  *r»  (%  +  O  d{ J 

A$’  *  cos6£J n,(tp+{)t flf)efdt+s*" 6<J *»•(*,+ i)efl9)d€ 


Aff  »  cos  6, 


117 


TABLE  H  -  Continued 


m  (ep+f)df +ep(eti+ep)J  md£  -  cos  djYtf/r  >*•«,) J  m(e0  +  e)d£ 

ir>6ct t,+eA  J  me  6 f  U  t 


+  Jt/U 


3{t  *  (e?  +tD)  J  m  (ep  +  4)d( 

CC  *  *']  ( *•(*,*■  VeU  +  *0*6^/  m(tr+4)e6g4{  +*1*6^  m(cp+4)(e0+e)d( 

*  *•  *4  Ip 

ca  mJ  ”>(*r+*'>(e,+S)dt 

A¥0^  *  -sin&^J  m  (e#  *■  f)(e0  *■  e)d(  -  cos  0* J , '/r?(e#-rf)e0fd4 

A 00  w  J mt0(to+Ze)d4  + J /ud4+m(y*+J*)  *■ J /n^vrftjj  +  jw^jj  Jdf 
Aj£  •J  *» [e0e(ic?K*+  (c.+e)/"']  dl 
A ef  m  fo  mec0Sjf)d4  +  J*  /ife*dt 

B9i  *  sin 6^(e0~ep)J (e9  / e)rr>d£  +  cos  Ot(eo-e0)J me  0,d( 

+  COS  &t  St  'n  d>'  u  Ke*md4  +  zJ*me0edl  +  Fr>  (tf.  j  *)] 

Z&*  \J meeo&i<Li  + J  Y’Stn  Z6fdB  ~  y  j| 


4- COS  i 


" a 

re 


t-sinze^J  ~C0SZ6(d{ 


118 


*ijwr 


TABLE  IE  -  Continued 


B6pt  -  Jm(e#tf)(?0  -t')d(  +cos9ijfn(ttt+i)(e,i-*)dt  -sme^Jm(eH*i)t&i  d{ 


CjZ}  -  26.  *<  +]*»•*£ dt  -  ifmrtgkffrPdt  j 

-Zsin  t6±  |  tff-  d  1 1  *  ^  coj  '.■[tjUlw-w]* 

*■  2jr*efl',k,',di  | 

^  -  z  fMi/<[pjrcw]  -  ^.j^,  c’Kr,Cd(\y'r> 

£J'[  -  cos  26; J  me,  [(tTK'T-  (f!f"):\  *t 

+  sin  6;  cos  ejm  [(tf&WV]  d{ 

C$.  -  cos‘6t^J mt9C  *  J^-  f^d.  f| 

-  Jt/?  ^  cos  0t  [  f  me  e0  6^  f^'df  +  zJpL0ff^f)d{^ 

’  cosetf\(e«  +{)  [CC*  ~  ]  <*■* 

s£»  &■  J  *m  (cH  *  ()  |VjTV  C'f?8t  ]  di 

m  coS  °*Jom(*r  r  #  \Cr>&f  ~  CV.T  ]  d  / 

+  Sin8;J  m(e^  4-f)  [C+  j  d{ 

CC"  / 'J'rjy+wwy'Xdi 

j 


119 


TABLE  TV  -  Continued 


Jmdf  (CK'l'-iC)*}*?} 


cy  - 

c£*  =  fU‘Vr>+ 1'"’;'"’]  A(*W 

♦  -jUc'crdt] 

^■cos&^tn  Gl  Jm  / J! ^ ^d {+2 Jm9f  | (k.tS>fr> )*-(£*)*  A 

-  sin% \tj”6t  C'K'C'd*  +J”  (C)  *dS  J 

cJ7.  *  f**>(e*+f)et'<rU( 

*  */o 

*■  cos*9i  j Jfn(eteo)0ffyd4  " +  k'r  d-% J 

+  cos0is£n0i  f  fm(ete0)/J,,r'd4  *■ J'n(e.+2e) 0f  ^4^'dK - //ra^VjTW 
U0  Jo  Jo 

r* 

-  sinz&t  / me  d.  f^d,  £ 

Jo 

CT  -  J  »(e«+f)e/%<9W 

4  0 o 

+  co$20t  I frte&f  f  ~ Jme Z  j 

+  cos  6tSin  0t  j  J mc  C’CW'J  {  "  ^Jme  | 

-  *>,'*•  { jmcef  C%%vdt * Jle  CWU  i  j 


*  cos 4 


+  cosi 


+  cost 


F(r.r) 

ur 


/'/«[, ;f..,[cn 


120 


TABLE  H  -  Continued 


Fj%  *  -  VC;  J  J ifdi  +  f  A(4)eA  i 


where 


&(Z)  -  /  (%+Z)md( 

Jf 

^  <r.p  _  d&l  C  (r>  v  r 

<W  -J0  EB*  jf  *»  -zj-dz 

-  -cosO^J  rr?(Cn  +  f )  £  &g  dg  -  *<>”  9;  J  mie#  +  Z)  Cf+  cLi 

a??-  J^WYdt 

Cjr  •  Cos%  j  J 4  J^V’di] 

+  Sin  &t  cos  Q^Jmtce  6f  (  +  V Jp 6 t  j  -SinO-  f  j pf^-f^aU  | 

C$.  *  cos\  [  Jmte0  Gf  t{a)di  *  f 

+  cos  dL  Sin  e-^Jmc.e0  f^all  +■  J m-  Sj'cSfJ  -  si»  *&<  Jfi  &f  f^cL  $ 


FP, 


'*  t  j/r  dfi*] 


- *—  d4 

d<  d( 


v  ■  Alh-% 


121 


APPENDIX  IV 


A  DEFINITION  OF  THE  GENERALIZED  FORCES  RESULTING 
FROM  AERODYNAMIC  LOADINGS  AND  OTHER  RESTRAINTS 


In  a  fully  articulated  rotor,  the  external  forces  on  the  blades 
arise  from  (1)  aerodynamic  loadings,  (2)  the  drag  link  and  damper,  and 
(3)  the  pitch-control  link.  These  external  forces  must  be  expressed  as 
generalized  forces  or  moments  applied  to  the  generalized  degrees  of 
freedom.  It  is  the  purpose  of  this  appendix  to  define  the  aerodynamic 
transformations  used  in  this  study  and  tc  present  the  expressions  used 
to  compute  the  generalized  forces  designated  by  the  symbol  "GF". 

If  experimental  lift,  drag,  and  moment  data  are  used  to  determine 
the  aerodynamic  forces  on  a  spanwise  segment  of  the  rotor  blade,  the 
segment  angle  of  attack  and  Mach  number  are  evaluated  as  functions  of 
the  kinematics  of  the  rotor  and  the  velocity  induced  by  the  wake.  Existing 
experimental  data  were  curve-fitted  by  expressions  similar  to  those 
developed  by  Blankenship  and  Harvey  [Reference  1].  (The  expression 
yielding  the  offset  of  the  center  of  pressure  from  the  one-quarter  chord 
has  already  been  given  in  Appendix  II. )  By  specifying  the  angle  of  attack, 
j  ,  to  range  between  plus  180  degrees  and  minus  180  degrees,  the  lift 
and  drag  coefficients  were  evaluated  from  the  following  expressions: 
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and  the  airfoil  section  properties,  CtJLi4mO),  SimM(/JtO),S0 >&i,  and  are 
experimental  data  applicable  to  the  airfoil  section  employed  in  a  given  rotor. 

Since  lift  and  drag  are  defined  as  forces  perpendicular  and  par¬ 
allel,  respectively,  to  the  local  total  velocity  vector,  the  normal  and  chord- 
wise  forces  on  the  blade  are  given  by 

JT 
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where 


L j  "  P  bi 

"  **  uw  bj  * 

If,  however,  the  aerodynamic  loadings  are  determined  by  means  of  the 
aerodynamic  model  described  in  Appendix  II,  the  normal-force  loading, 
r  •  •  ,  and  the  chordwise-force  loading,  ,5,  ■  •  ,  are  obtained  directly 
from  the  analysis. 

In  order  to  obtain  (1)  the  drag  force  on  the  blade  parallel  to  the 
shaft  plane  and  (2)  the  pitching  moment  about  the  feathering  axis,  it  was 
convenient  to  resolve  the  normal  and  chordwise  forces  into  other  axis 
systems  and  to  define  the  moment  arms  to  the  feathering  axis  in  terms 
of  coordinates  based  on  an  axis  system  whose  </  axis  is  located  along 
the  chord  of  the  blade  root.  Figure  9  (see  page  23)  shows  that  "lift  and 
drag"  forces  defined  with  respect  to  the  y-jand  y„  axes,  respectively,  are; 


A 

m 

m 

A  a> 

Stn  An,  i,  j  An,  »,  j 

A 

A 

m 

A  » 

m 

stn 

COa 

- 

cos  &  ; 

m 

Stn 

- 

r 

Si*  ;  : 

COS  Qm  ■  ;  Dm  ■ 

where 


*, J  '  A'  +  *,j 

~  Q*.*  + 


By  utilizing  the  kinematic  derivation  procedure  outlined  in 
Appendix  I,  the  following  expressions  are  obtained  for  the  distance  of 
the  one-quarte.'  chord  point  from  the  feathering  axis 

+  I'M j  A  +  $*•,*,  j 

Ui.j  ’  e  [«.;  '•  *5,.;.;  ♦  1w]  ♦  O'm.i.i  *  r»j 

where  we  have  assumed  that  cos  6:  ■  f  and  sin  Q  ■  &; ■ 
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The  generalized  forces  and  moments  presented  below  were 
derived  by  using  a  procedure  similar  to  that  outlined  in  Appendix  V  of 
Reference  [8].  It  should  be  noted  that  the  concentrated  load  arising 
from  the  drag  link  was  assumed  (in  this  study)  to  act  on  that  portion  of 
the  blade  (inboard)  which  is  taken  to  be  infinitely  stiff  or  rigid.  Accord¬ 
ingly,  the  drag-damper  force  does  not  enter  into  the  generalized  forces 
applicable  to  the  blade-bending  degrees  of  freedom.  The  drag  damper 
does,  however,  contribute  to  the  generalized  moment  *.bout  the  lead-lag 
hinge. 


On  neglecting  the  effects  of  blade  bending,  lead-lag  displacement, 
and  flapping  displacement  on  the  moment  arm  to  the  rotor  shaft,  we  have 
that  the  aerodynamic  torque,  £>4  ,  produced  by  blade,  /  ,  is 

6,  i  -  -  /  (e„  *  f)  D„.i  U)d(  . 

The  total  external  torque  (or  generalized  force)  about  the  rotor  shaft  is 

*/s 

QM  +Z  Qa  . 

where  Q(t)  is  the  mechanical  torque  produced  by  the  helicopter  power 
plant. 


On  neglecting  the  effects  of  blade  bending  and  lead -lag  displace¬ 
ment  on  the  moment  arm  to  the  flap  hinge,  the  generalized  flapping  moment 
can  be  expressed  as 

uo 

Similarly,  the  generalized  lead-lag  moment  can  be  expressed  as 

«■ C  -  Ki Nr  C... 

where 

is  the  Urag-damper  moment  per  unit  lead-lag  velocity, 
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The  generalized  moment  about  the  feathering  axis  can  te  expressed 
as 

mJo  (VL.i (O +  #,[%(*)  -OmJ] 

where  Mg  *']  *■  the  pitching  moment  caused  by  coi.ipreBsion  or 

tension  in  the  pitch  link  and  £(-ft)is  the  effective  swashplate  angle  for 
blade  i ,  as  caused  by  the  combined  collective-  and  cyclic -pitch  settings. 

It  should  be  noted  that  the  pitching  moment  can  be  defined  to  include  other 
torques  such  as  would  be  caused  by  hub  properties  that  yield  pitching 
moments  due  to  displacements  about  the  flapping  hinge. 

The  generalized  force  for  the  rth  mode  of  coupled  flatwise  and 
chordwise  bending  is 

and  the  generalized  torque  for  the  J  **  torsion  mode  is 

where  ^  is  the  distance  between  the  elastic  axis  and  the  one -quarter  chord 
and  is  the  aerodynamic  pitching  moment  about  the  one -quarter  chord. 
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